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SECTION  I 


INTRODUCTION 


A.  OBJECTIVE 

The  objectives  of  this  research  are  to  determine; 

1)  scaling  laws  for  explosion-induced  apparent  crater 
dimensions  of  volume,  radius,  and  depth  for  dry  and 
partially  saturated  sand  for  buried,  fully  coupled 
explosives ; 

2)  the  effect  of  the  depth  of  charge  burial  (DOB)  on 
crater  formation; 

3)  to  what  extent  desaturation  occurs  in  beach  sand 
during  centrifugation; 

4)  the  effect  of  moisture  on  the  apparent  crater 

dimensions ; 

5)  the  effect  of  moisture  on  the  apparent  crater 

aspect  ratio; 

6)  and  the  effect  cf  the  Coriolis  acceleration  on  the 
shape  of  the  apparent  crater. 

For  the  first  objective,  the  apparent  crater  dimensions 
are  analyzed  as  dimensional  terms  and  as  dimensionless  pi 
terms.  In  the  analy  is  of  dimensional  terms,  the  apparent 


crater  dimensions  in  dry  and  partiai3.y  saturated  sand  are 
plotted  against  the  TNT  equivalent  weight  of  the  explosive 
mass.  The  slope  of  the  best-fit  lines,  the  yield  exponents, 
are  compared  to  determine  the  similarity  between  scaling  laws. 
Comparisons  of  the  yield  exponents  derived  in  this  subscale 
centrifuge  research  are  also  made  with  the  yield  exponents 
derived  from:  1)  other  centrifuge  tests  (Goodings  et  al., 
1988);  Schmidt  and  Holsapple,  1978  and  1980);  2)  field 
explosive  tests  by  the  U.S.  Army  Engineer  Waterways  Experiment 
Station  (Strange  et  al.,  1961;  Roolce  et  al . ,  1974);  and  3) 
crater  scaling  relations  given  in  the  Air  Force  Protective 
Construction  Design  Manual  and  the  Army  de.sign  manual. 
Fundamentals  of  Protective  Design  for  Conventional  Weapons. 

The  dimensionless  pi  terms  used  to  analyze  crater 
dimensions  produced  by  buried  charges  are  those  derived  by 
Schmidt  and  Holsapple  (1977)  .  In  the  analysis  of  the  pi 
terms,  a  direct  comparison  is  made  between  the  consistency  of 
the  apparent  crater  dimensions  with  increasing  charge  size  in 
dry  sand  and  in  partially  saturated  sand. 

B.  BACKGROUND 

Theoretical  research  in  cratering  consists  of  numerical 
methods  that  either  combine  a  cavity  growth  and  surface 
mounding  model  with  a  free-fall  throw  out  model,  or  use  quasi¬ 
static  calculations  to  treat  cratering  as  an  earth  pressure 
problem  (Dillon,  1972).  The  experimental  research  has 
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concentrated  on  the  development  of  scaling  relations  that 
allow  the  crater  dimensions  resulting  from  one  charge  size  to 
be  scaled  to  another.  The  vast  majority  of  the  cratering  data 
generated  through  experimental  research  comes  from  full-scale 
field  tests.  Unfortunately,  much  of  the  field  test  data 
includes  little,  if  any,  information  on  the  soil  properties. 

C.  SIGNIFICANCE 

Cohesionless  soils  are  commonly  placed  and  compacted 
around  and  over  strategic  structures.  The  influence  of  the 
moisture  content  of  cohesionless  soils  on  ground  shock  and 
cratering  from  conventional  weapons  detonations  is  unknown. 
Unsccuratecl  soil  mechanics  theory  predicts  saturation  level 
will  influence  soil  behavior  (Fredlund,  1986) . 

Split-Hopkinson  Pressure  Bar  tests  conducted  by  Ross  et 
al .  (1989)  and  Pierce  (1989)  indicate  that,  in  sands,  the 
saturation  during  compaction  influences  stress  transmission. 
Uniaxial  explosive  tests  conducted  on  the  centrifuge  by 
Charlie  and  Walsh  (1990)  indicate  trends  similar  to  those  of 
the  Split-Hopkinson  Pressure  Bar- 

Centrifuge  test  results  on  the  material  strength  model 
proposed  by  Schmidt  and  Holsapple  (1978,  1980)  indicate  that 
explosive  induced  crater  dimensions  scale  differently  for 
soils  with  cohesion  (moist  alluvium  and  clays)  than  for  dry 
frictional  cohesionless  soils  (sands  and  gravels)  .  The 
difference  in  scaling  occurs  when  cohesion  is  large  compared 
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to  the  total  shear  strength.  Therefore,  the  test  results  to 
date  indicate  that  the  moisture  content  of  compacted  sands  may 
influence  the  conventional  weapon  effects  on  airbase 
structures . 

D.  APPROACH 

The  observational  method  of  analysis  through  centrifuge 
modeling  is  used  to  study  explosion-induced  cratering. 
Centrifuge  modeling  allows  soil  properties  to  be  well- 
controlled  and  provides  the  gravitational  field  required  for 
similarity  between  model  and  prototype.  Stiffness  and 
strength  of  dry  cohesionless  soils  are  derived  from  the  self¬ 
weight  stresses.  Therefore,  only  by  increasing  the 
gravitational  acceleration  can  the  self-weight  stresses  be 
equivalent  in  the  model  and  prototype. 

The  scale  of  the  model  determines  the  gravitational 
acce?..eration  required  to  equate  the  soil  properties  in  the 
model  and  prototype.  A  1/N  scale  model  will  require  N  times 
the  earth's  gravitation  il  constant  (G  =  9.81  m/s^)  to  maintain 
similarity  between  the  model  and  prototype.  For  similarity  of 
explosive  energy,  a  1/N  scale  model  at  N  times  the  earth's 
gravitational  constant  will  require  1/N^  times  the  explosive 
mass  of  the  prototype.  In  a  1/N  scale  crater,  the  depth  and 
diameter  scale  with  1/N,  whereas  the  crater  volume  scales  with 
1/N\  Centrifuge  tests  will  be  conducted  at  four  G-levels  to 
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verify  the  application  of  these  and  other  scaling  laws  to 
unsaturated  sands. 
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SECTION  II 


FULL  SCALE  EXPLOSION-INDUCED  CRATERING 

A.  CRATER  DIMENSIONS 

Figure  2.1  shows  the  dimensions  used  to  describe  a 
crater.  The  ejected  material  deposited  on  the  ground  surface 
surrounding  the  crater  by  the  blast  is  referred  to  as  the 
ejecta.  The  airborne  material  that  falls  back  into  a  newly 
formed  crater  is  referred  to  as  fallback.  The  subscripts  a 
and  t  indicate  the  apparent  and  true  crater  dimensions, 
respectively.  The  apparent  crater  is  the  visible  crater  below 
the  level  of  the  original  ground  surface.  The  true  crater  is 
defined  by  the  boundary  below  the  original  ground  surface 
between  the  fallback  and  the  crushed,  fractured  material  that 
did  not  experience  significant  displacement  or  disarrangement 
(Dillon,  1972)  . 

The  crushed  and  fractured  material  that  does  not  undergo 
significant  displacement  or  disarrangement  makes  up  the 
rupture  zone.  The  plastic  zone  defines  the  region  in  which 
penrianent  deformation  has  occurred  without  crushing  or 
fracturing.  The  deformation  resembles  a  flow  of  the  material. 
Within  the  plastic  zone  is  a  shear  boundary  beyond  which  no 
shearing  of  the  material  is  assumed  to  occur.  The  elastic 


r^=  True  Crater  Radius 
d^=  True  Crater  Depth 


Figure  2.1  Half-Crater  Profile  taken  about  a  Vertical 
Centerline. 
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zone  contains  material  that  has  not  undergone  permanent 
deformation  (Rooke  et  al.,  1974). 

B.  CRATERING  TERMINOLOGY 

Centrifuge  crater  dimensions  are  described  as  measured, 
scaled,  or  normalized.  Scaled  crater  dimensions  are  measured 
model  dimensions  multiplied  by  N,  where  N  is  a  multiple  of  the 
earth's  gravitational  acceleration  (G  =  9.81  m/s^).  Thus, 
scaled  crater  dimensions  are  the  prototype  crater  dimensions. 

In  crater  modeling,  dimensions  are  divided  by  w'’  -o 
eliminate  the  dependence  on  explosWe  energy.  Charge  mass  (W) 
is  commonly  used  to  represent  explosive  energy  because  for  a 
given  chemical  e.xplosive,  charge  mass  is  proportional  to 
explosive  energy.  With  all  other  variables  held  constant,  the 
charge  mass  exponent,  n,  is  referred  to  as  the  yield  exponent. 
When  the  yield  exponent  is  unknown,  a  value  of  1/3  has  been 
found  to  be  a  good  approximation  for  dimensions  of  radius  and 
depth  (Strange  et  al . ,  1961)  .  The  procedure  of  dividing  by  W" 
is  referred  to  as  normalization.  Thus,  the  normalized  crater 
dimensions  are  scaled  dimensions  divided  by  the  cube  root  of 
the  charge  mass. 

Depth  of  burial  (DOB)  significantly  effects  crater  size. 
To  determine  the  functional  relationship  between  DOB  and 
crater  size,  normalized  crater  dimensions  are  plotted  against 
normalized  DOB  (K) •  As  with  crater  radius  and  depth,  DOB  is 
divided  by  W''^.  Py  convention,  a  positive  DOB  refers  to  a 
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buried  charge,  and  a  negative  DOB  refers  to  an  airblast.  The 
height  of  burst  (HOB)  is  also  used  to  denote  the  distance  of 
an  airblast  from  the  ground  surface. 

C.  CRATERING  MECHANISMS 

Near  surface  blasts  produce  craters.  The  size  and  shape 
of  the  crater  are  controlled  by  the  energy  of  the  explosion, 
the  properties  of  the  soil  or  rock,  DOB,  and  gravity  (Dillon, 
1972) .  For  all  other  variables  held  constant,  the  greater  the 
energy  release,  the  larger  the  crater.  Conversely,  the 
greater  the  shear  strength  of  the  soil,  the  smaller  the  crater 
for  a  given  release  of  energy  and  DOB. 

The  fraction  of  the  released  energy  that  impacts  the 
medium  is  a  function  of  DOB.  Initially  as  DOB  increases,  the 
crater  becomes  larger  because  more  energy  is  directed  into  the 
medium.  Then  beyond  a  certain  DOB,  the  impact  of  the  energy 
released  on  the  surface  becomes  less  and  less  and  the  size  of 
the  crater  begins  to  decrease.  The  optimum  DOB  is  the  DOB 
that  maximizes  a  crater  dimension.  The  optimum  DOB  for  crater 
volume,  radius,  or  depth  are  not  necessarily  the  same  (Dillon, 
1972)  . 

Gravity  is  a  significant  factor  in  crater  formation. 
First,  the  explosive  energy  must  overcome  the  gravitational 
force  to  excavate  the  soil  or  rock.  Second,  gravity  controls 
the  distribution  of  fallback.  Third,  gravity  controls  the 
properties  of  frictional  materials 
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There  are  three  mechanisms  for  nonnuclear  crater 
production.  In  chronological  order  these  mechanisms  are  the 
formation  of  a  shock  front,  explosive  gas  acceleration  and 
gravity.  The  shock  front  is  the  pressure  discontinuity 
developed  from  the  explosive  detonation.  The  shock  front 
impacts  the  soil  first  and  causes  crushing,  compaction, 
fracturing,  and  plastic  deformation  of  the  soil  surrounding 
the  blast  cavity.  Compression  waves  are  reflected  at  the 
ground  surface  as  a  tensile  wave.  When  the  magnitude  of  the 
reflected  wave  exceeds  the  tensile  strength  of  rhe  material, 
spalling  occurs  (Dillon,  1972)  . 

During  the  gas-acceleration  phase,  the  crushed  and 
fractured  material  is  accelerated  by  the  adiabatic  expansion 
of  the  trapped  gases.  The  magnitude  of  the  acceleration 
vector  determines  the  flight  of  the  material  particles. 
Following  the  gas  expansion,  the  gravitational  force  pulls  the 
airborne  particles  down  into  and  around  the  crater. 

D.  HISTORICAL  CRATERING  INVESTIGATION  AT  1  G 

Initial  research  on  crater  prediction  and  scaling  relied 
on  field  test  data.  The  goal  was  to  derive  the  functional 
relationship  between  apparent  crater  radius  (r^),  depth  (d^) 
and  explosive  energy  (W)  .  These  functional  relationships  are 
commonly  referred  to  as  scaling  laws,  and  are  named  after  the 
yield  exponent.  For  example,  the  cube  root  scaling  law 
describes  a  condition  in  which  the  linear  dimensions  of  the 
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model  crater  scale  to  prototype  dimensions  with  the  cube  root 
of  the  explosive  energy  (r,  d  «  .  in  cube-root  scaling, 

crater  volume  is  proportional  to  explosive  energy  (V  «  W) 
because  volume  is  proportional  to  the  cube  of  a  linear 
dimension . 

Other  variables  included  in  the  research  on  cratering 
were  DOB,  soil  strength  (a)  expressed  in  dimensions  of  stress, 
and  the  total  unit  weight  of  the  soil  (K) .  The  soil  strength 
variable,  d,  was  used  to  represent  the  constitutive  effects  of 
the  medium  on  the  crater  dimension,  and  the  unit  weight 
variable,  K,  was  used  to  represent  the  effect  of  gravity. 
Through  dimensional  analysis  (a  variable  manipulation 
technique  used  to  generate  the  conditions  for  model  and 
prototype  similarity  and  discussed  in  section  IV)  the 
functional  relationship  including  all  of  the  above  named 
variables  was  found  to  be 


dob  o-''^DOB  ’  K'-^^DOB 


(2.1) 


The  constitutive  effects  term  (W‘'^/a‘'^DOB)  follows  the  cube 
root  scaling  law,  and  the  gravitational  effects  term 
{W-'VK‘''‘D0B)  follows  the  quarter-root  scaling  law  (Baker  et 
al . ,  1973)  . 

To  reduce  the  three-dimensional  functional  relationship 
to  two  dimensions,  several  researchers  ignored  either  the 


constitutive  effects  or  the  gravitational  effects.  Chabai 
<1965)  recognized  that  neither  the  constitutive  effects  term 
nor  the  gravitational  effects  term  alone  were  sufficient  to 
describe  apparent  crater  radius.  Chabai  (1965)  stated  that, 
for  small  differences  in  explosive  energy,  the  cube-root 
scaling  law  applied,  but  that  for  large  differences  in 
explosive  energy  the  quarter-root  scaling  law  provided  a 
better  estimate  of  apparent  crater  dimensions.  Earlier,  in  a 
regression  analysis,  Chabai  (1959)  found  the  r/DOB  term  could 
be  approximated  by  the  function  W*'^  VdOB.  The  fact  that  the 
exponent  for  the  charge  mass  was  between  1/3  and  1/4  was  ta)<en 
by  Baker  et  al .  (1981)  as  an  indication  that  in  cratering, 
both  the  constitutive  effects  and  the  gravitational  effects 
are  significant. 

In  1960  and  1961,  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  published  a  two-volume  technical 
report  entitled  "Cratering  from  High  Explosive  Charges."  The 
report  shows  cube-root  scaling  best  describes  crater  radius 
and  a  yield  exponent  of  0.3  best  describes  crater  depth 
(Strange  et  al .  ,  1961)  .  However,  the  authors  of  the  WES 
report  state  that  the  cube-root  scaling  law  is  only  accurate 
for  charges  ranging  from  0.5  kg  to  450  kg  of  TNT.  They 
attribute  the  deviations  from  the  principles  of  similitude  to 
the  dependence  of  soil  strength  on  gravity  and  the  difficulty 
in  scaling  soil  density,  particle  size,  and  soil  strength 
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properties  in  proportion  to  the  ratio  of  the  cube  root  of  the 
model  and  prototype  charge  masses. 

The  data  are  also  normalized  to  allow  comparison  of  the 
various  test  series  and  to  form  prediction  curves .  The 
prediction  curves  show  normalized  crater  radius  (r/W''^^)  and 
normalized  crater  depth  (d/W^'^^)  plotted  against  normalized  DOB 
(X=DOB/W^^^)  .  The  maximum  apparent  crater  radius  and  depth  are 
shown  to  occur  when  1.0  <  X  <  1.5,  and  the  maximum  true  crater 
radius  and  depth  are  shown  to  occur  near  X  =  2.0. 

Rooke  et  al .  (1974)  updated  the  WES  report.  Data  from 
nuclear  explosions  are  included  with  the  high  explosive  data. 
The  data  were  also  separated  by  moisture  content,  as  well  au 
by  medium  and  DOB. 

Dillon  (1972)  performed  a  regression  analysis  on  existing 
crater  data  to  develop  prediction  equation^  for  apparent 
crater  dimensions  as  a  function  of  DOB,  and  to  determine  the 
effect  of  material  properties  on  the  apparent  crater 
dimensions.  Dillon  (1972)  used  10  nuclear  events  with  yields 
ranging  from  308,000  kg  to  91,000,000  ky  TNT  equivalent,  and 
200  high  explosive  events  with  yields  ranging  from  0.45  kg  to 
454,000  kg  of  TNT.  The  densities  of  the  materials  ranged  from 
960  kg/m^  to  2,720  kg/m-*.  A  bell  shaped  curve  was  selected  for 
the  regression  model  and  a  value  of  5/) 6,  the  average  of  1/3 
and  7/24,  was  used  as  the  yield  exponent. 

Dillon  (1972)  states  the  "best  regression  equations" 
occurred  when  the  soil  property  variables  consisted  of  total 
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unit  weight,  percent  saturation,  specific  gravity,  shearing 
resistance,  and  seismic  velocity.  Of  these  variables,  percent 
saturation  and  total  unit  weight  appeared  to  be  the  most 
significant.  At  the  "optimum"  (undefined,  presumably  for 
crater  diameter)  DOB,  "larger"  (undefined,  presumably  also  in 
terms  of  crater  diameter)  craters  were  produced  at  low  and 
high  levels  of  saturation  in  granular  soils.  The  "largest." 
craters  at  the  "optimum"  DOB  occurred  in  dry  soil . 


SECTION  III 


CENTRIFUGE  MODELING 

A.  PRINCIPLE  OF  CENTRIFUGE  OPERATION 

Centrifuge  modeling  was  first  applied  to  geomechanical 
testing  in  1931  by  P .  B.  Bucky  at  Columbia  University,  and  in 
1932  by  G.  1.  Pokrovsky  and  N.  N.  Davidenkov  in  the  U.S.S.R. 
The  basic  components  of  a  centrifuge  used  for  geomechanical 
testing  are  shown  in  Figure  3.1a.  The  centrifuge  increa.ses 
the  acceleration  applied  to  a  specimen  on  the  payload  platform 
by  rotating  in  the  horizontal  plane.  Centripetal 

acceleration,  a^,,  is  directed  radially  toward  the  center  of 
rotation  and  is  defined  as 


where  O)  is  the  angular  velocity  in  radians  per  second  and  r  is 
the  radial  distance  to  the  specimen.  The  centripetal  force, 
Fc  is 

(3.2) 

where  M  is  the  specimen  mass.  Centripetal  force  is  present 
when  the  frame  of  reference  is  outside  the  rotating  frame  and 
is  therefore  defined  as  a  real  force.  The  equal  and 


Rotor  Arms 


Figure  3.1  Centrifuge  Schematic:  a)  Centrifuge  Components; 

b)  Payload  Platform  Rotation  During  Centrifuge 
Operation  and  Force  Direction-s. 
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oppositely  directed  force  is  the  centrifugal  force  (Figure 
3.1b).  The  centrifugal  force  acts  only  in  the  rotating  frame 
and  is  therefore  a  pseudo  force. 


B.  ASSUMPTIONS  IN  CENTRIFUGE  MODELING 

The  centrifuge  modeling  of  a  process  or  event  assumes:  1) 
the  model  is  an  accurate  scaled  version  of  the  prototype;  2) 
a  1/N  scale  model  at  N  G  behaves  like  the  prototype  at  1  G; 
and  3)  the  centrifuge  produces  an  N  G  gravitational  field  that 
is  constant  throughout  the  specimen  (Joseph  et  al.,  1988) . 

The  accuracy  of  the  model  depends  on  the  scaling 
relations  between  the  model  and  the  prototype.  A  scaling 
relation  is  the  expression  relating  the  prototype  variable  at 
1  G  to  the  model  variable  at  N  G.  For  example,  the  scaling 
relation  for  a  linear  dimension,  L,  is 

=  (3.3) 

where  the  subscripts  p  and  m  refer  to  prototype  and  model,  and 
the  constant  N  is  the  G-level .  Table  3.1  shows  the  scaling 
relationships  used  in  centrifuge  modeling. 

The  second  assumption  in  centrifuge  modeling  states  that 
whatever  happens  in  the  prototype  will  happen  in  the  model. 
Dynamic  similarity  between  model  and  prototype  exists  only 
when  the  nondimensional  parameters  or  pi  terms  formed  from 
prototype  variables  equal  the  equivalent  nondimensional 
parameters  formed  from  model  variables.  For  example,  when 
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TABLE  3.1.  SCALING  RELATIONS  FOR  CENTRIFUGE  MODELING. 


Full 

Scale 

Centrifuge 

Quantitv  (Prot 

:otvDe) 

Model  at  N  G 

Linear  Dimension 

1 

1/N 

Area 

1 

1/N^ 

Volume 

1 

1/N^ 

Time 

In  Dynamic  Terms 

1 

1/N 

In  Diffusion  Terms 

1 

1/N^ 

In  Viscous  Flow  Cases 

1 

1 

Velocity  (Distance/Time) 

1 

1 

Acceleration  (Di  stance/Time^) 

1 

N 

Mass 

1 

l/N^ 

Force 

1 

1/N^ 

Energy 

1 

1/N^ 

Stress  (Force/Area) 

1 

1 

Strain  (Displacement /Uni t  Length) 

1 

1 

Density 

1 

1 

Energy  Density 

1 

1 

Frequency 

In  Dynamic  Problems 

1 

N 

n  =  =  7i„  (3.4) 


where  V  is  apparent  crater  volume,  p  is  soil  density,  and  W  is 
charge  mass.  Therefore,  the  development  of  a  good  model 
depends  on  the  identification  of  the  controlling  variables  and 
the  generation  of  the  associated  pi  terms. 

There  are  three  methods  used  to  generate  the  pi  terms  for 
a  physical  phenomenon  such  as  an  explosion.  If  the  physical 
factors  contro] ling  the  phenomenon  can  be  identified,  then  a 


governing  differential  equation  can  often  be  written  to 
describe  the  phenomenon.  The  use  of  a  governing  equation  is 
the  first  method  that  can  be  used  to  insure  the  existence  of 
similarity . 

Pi  terms  are  created  from  a  governing  equation  through  a 
normalization  procedure.  Kline  (1980)  defines  noirmalization 
as  "making  the  governing  equations  and  conditions 
nondimensional  in  terms  of  nondimensional  variables  of 
standard  magnitude."  The  procedural  steps  are  first  to  make 
the  variables  dimensionless  (divide  each  variable  by  itself 
defined  at  the  limit  of  the  domain  of  interest) ,  substitute 
the  dimensionless  variables  into  the  governing  equation^  and 
lastly  to  make  the  coefficient  of  each  term  dimensionless 
(divide  each  coefficient  by  the  coefficient  of  one  of  the 
terms).  The  remaining  coefficients  are  the  pi  terms. 

The  first  step  is  the  more  difficult  because  the 
appropriate  scales  of  the  problem  must  be  known.  The  main 
advantage  of  using  a  governing  equation  is  that  all  the 
controlling  variables  are  known  to  be  included  in  the 
analysis.  The  list  of  pi  terms  is  complete  and  physically 
meaningful . 

The  second  method  used  to  create  the  pi  terms  for  a 
physical  phenomenon  is  the  method  of  similitude.  The  method 
of  similitude  involves  identifying  the  forces  controlling  the 
phenomenon,  defining  the  forces  in  the  term  of  relevant 
variables,  and  creating  force  ratios.  The  force  ratios  are 
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the  pi  terms.  Thus,  the  method  of  similitude  is  useful  when 
the  physical  process  is  not  well  understood  but  the 
controlling  forces  are. 

When  the  physical  process  is  not  well-understood  and  the 
forces  cannot  be  well-defined,  the  method  of  dimensional 
analysis  is  used  to  generate  a  list  of  pi  terms.  Dimensional 
analysis  is  an  extension  of  the  Buckingham  pi  theorem 
(Buckingham,  1914)  which  states  that  given  a  relation  among  m 
parameters  of  the  form 

Pi,  Pj.  .  .  .  .  p„  )  =  0  (3.5) 

an  equivalent  relation  among  n  terms  of  dimensionless 
parameters,  pi  terms,  can  be  found  of  the  form 

F2i  Ttz'  •  •  •  '  ^  °  (3.6) 

where  n  is  usually  m  less  the  number  of  primary  units 
represented  by  the  parameter. 

The  dimensional  analysis  of  a  physical  phenomenon  begins 
with  the  listing  of  ail  relevant  variables  expressed  in 
primary  units.  There  are  two  systems  of  primary  units;  force, 
length,  and  time  or  mass,  length,  and  time.  The  mass-length¬ 
time  system  should  be  selected  for  centrifuge  modeling  as 
force,  unlike  mass,  varies  with  G-level. 

Dimensionless  combinations  of  the  variables,  the  pi 
terms,  are  created  through  algebraic  cr  matrix  manipulation. 
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The  pi  terms  are  complete  for  the  variables  selected  as  long 
as  each  pi  term  is  independent  of  the  others,  and  any  pi  term 
not  generated  by  the  analysis  can  be  formed  by  taking  the 
product  of  two  or  more  of  the  terms  raised  to  some  power. 
However,  the  pi  terms  are  not  absolute  if  important  variables 
are  missing.  Thus,  the  possible  exclusion  of  important 
variables  is  the  main  weakness  of  dimensional  analysis. 

Thare  are  two  significant  advantages  of  using  pi  terms  to 
study  a  physical  phenomenon  other  than  to  insure  similarity. 
The  first  is  that  combinations  of  variables  allows  certain 
variables  to  vary  between  model  and  prototype.  Only  the 
combinations  themselves  must  remain  equal .  Second,  the 
process  of  generating  the  pi  terms  may  give  physical  meaning 
to  seemingly  unrelated  variables. 

The  third  assumption  of  centrifuge  modeling  is  that  the 
uniform  nature  of  the  earth' s  gravitational  field  at  the 
surface  can  be  accurately  reproduced  at  higher  G-levels. 


C.  SOURCES  OF  ERROR  IN  CENTRIFUGE  MODELING 

One  obvicus  source  of  error  in  centrifuge  modeling  is  the 
boundary  effect.  The  specimen  has  to  be  placed  in  a  container 
of  one  form  or  another.  However,  boundary  effects  can  usually 
be  shown  to  be  negligible  in  comparison  to  the  measured  event. 

A  second  source  of  error  occurs  when  one  or  more  of  the 
scaling  relations  can  not  be  met.  Under  these  conditions  the 
model  becomes  distorted.  Limited  distortion  of  a  model  is 


21 


acceptable  as  long  as  the  physical  phenomenon  is  not  affected, 
or  the  test  results  are  corrected  for  the  distortion. 

Distorted  models  occur  in  geotechnical  centrifuge 
modeling  when  the  prototype  sand  is  used  in  the  model.  No 
attempt  is  made  to  scale  the  particle  grain  size,  because  to 
do  so  would  necessitate  the  use  of  clay-sized  particles  having 
different  material  properties.  The  difficulty  with  not  using 
a  scaled  grain  size  in  the  model  is  that  the  linear  dimensions 
of  the  sand  grains  scale  to  a  larger  material  in  the 
prototype.  However,  as  long  as  the  grain  size  distortion  does 
not  affect  other  areas  of  the  model  such  as  the  number  of 
contact  points  with  a  structure  or  instrument,  the  model 
distortion  is  considered  to  be  of  secondary  importance.  The 
primary  consideration  is  that  the  behavior  of  the  soil  be  the 
same  in  the  model  and  prototype. 

Through  a  similarity  analysis  of  cratering  using  the 
thermomechanical  response  of  a  continuum,  Schmidt  and 
Holsapple  (1980)  showea  that  dynamic  similarity  between  tests 
in  the  same  material  wili  occur  as  long  as  the  constitutive 
behavior  of  the  material  is  independent  of  the  scale  factors 
for  time  and  size. 

Constitutive  behavior  i.idependent  of  the  scale  factors 
for  time  is  rate  independent.  Thus,  sand  used  in  centrifuge 
modeling  should  ideally  be  rate-independent.  However,  Felice 
et  al .  (1985)  showed  sand  exhibits  strain-rate  dependence  when 
the  strain  is  greater  than  the  initial  gas  porosity.  Farr 
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(1990)  measured  20  to  100  percent  increases  in  constrained 
modulus  due  to  differences  in  loading  rate.  In  addition, 
recent  Split-Hopkinson  Pressure  Bar  tests  show  a  strain-rate 
dependence  in  moist  Eglin  sand^.  Thus,  additional  research 
is  required  to  determine  to  what  extent  rate  dependence 
affects  the  conditions  of  similarity. 

A  third  source  of  error  lies  in  the  use  of  centripetal 
acceleration  to  increase  the  G-level .  Centripetal 
acceleration  varies  with  distance  from  the  axis  of  rotation. 
Thus,  within  a  cylindrical  specimen,  lines  of  constant 
acceleration  follow  lines  of  equal  radii  (Figure  3.2). 

The  linear  increase  in  G-level  with  specimen  depth 
creates  a  nonlinear  stress  distribution  in  the  model  (Figure 
3.3)  .  Whereas,  a  uniform  gravitational  field  in  the  prototype 
leads  to  self-weight  stresses  that  increase  linearly  with 
depth.  Therefore,  the  model  stresses  equal  the  prototype 
stresses  at  only  one  depth.  Above  this  depth,  the  model 
stresses  are  lower  than  the  prototype  stresses,  while  below 
this  depth,  the  model  stresses  are  higher  than  the  prototype 
stresses.  The  error  associated  with  this  stress  difference  is 
a  function  of  the  ratio  between  the  model  depth  and  the  radius 
of  the  centrifuge  rotor  arm  (Joseph  et  al.,  1988).  A  long 
rotor  arm  and  a  shallow  model  can  be  used  to  minimize  the 
error . 

^Personal  Communication,  C.  Allen  Ross,  Professor  Emeritus, 
University  of  Florida,  Gainesville,  Florida;  Consultant,  HQ 
AFCESA/RACS,  Tyndall  AFB,  Florida. 
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Lines  of  Constant 


Variation  of  Centripetal  Acceleration 
Within  a  Specimen,  Side  View.  (NTS) 


Figure  3.2 


Soil  Pressure 


Figure  3.3  Model  Stress  Error  with  Depth 
(Joseph  et  al . ,  1988  )  . 
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A  second  acceleration,  the  Coriolis  acceleration,  is 
added  to  the  system  when  the  blast  occurs.  The  Coriolis 
acceleration  is  defined  as 

=  2q  XV  (3.7) 

where  O)  is  angular  velocity,  and  v  is  particle  velocity. 
Therefore,  to  determine  the  movement  of  particles  in  a 
rotating  frame,  the  Coriolis  acceleration  must  be  included. 
In  the  centrifuge  environment,  the  change  in  particle  movement 
due  to  the  Coriolis  acceleration  is  referred  to  as  the 
Coriolis  effect. 

Lastly,  the  earth's  gravitational  acceleration  cannot  be 
overlooked  at  low  G-levels.  The  acceleration  on  the  specimen 
is  the  resultant  of  the  centripetal  and  gravitational 
acceleration : 


conctipotal 


^  gzavltacional 


(3.8) 


In  summary,  the  centrifuge  does  not  meet  exactly  all  of 
the  assumptions  made  to  justify  its  use.  However,  the  errors 
can  usually  be  made  insignificant.  The  container  size  should 
be  large  enough  to  reduce  boundary  effects.  The  particle  size 
should  be  small  enough  to  maintain  sufficient  points  of 
contact,  but  not  so  small  as  to  alter  the  material  properties 
of  the  soil.  The  model  size  should  be  limited  to  some 
fraction  of  the  rotor  arm  length,  and  tests  should  be  run  at 
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speeds  high  enough  to  make  the  gravitational  acceleration  a 
small  percentage  of  the  resultant. 

D.  MODELING  OF  MODELS 

The  modeling  of  models  technique  refers  to  the  comparison 
of  .similar  centrifuge  tests  at  G-levels  bounding  the  domain  of 
interest.  If  the  measurements  from  each  test  are 

approximately  equal,  then  the  scaling  relations  are  shown  to 
be  accurate  over  the  domain.  No  scale  effects,  errors 
introduced  through  the  scaling  process,  are  present  and  the 
model  should  accurately  simulate  the  prototype.  The  advantage 
of  a  test  series  that  includes  the  modeling  of  models  is  that 
centrifuge  modeling  can  be  performed  without  a  prototype. 
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SECTION  IV 


CRATERING  ON  THE  CENTRIFUGE 

A,  DIMENSIONAL  ANALYSIS  -  PI  TERMS 
1 .  Surface  Detonations 

Schmidt  and  Holsapple  (1978)  performed  a  series  of  tests 
in  dry  Ottawa  sand  to  verify  the  similarity  relations 

generated  by  their  thermomechanical  analysis  of  a  continuum, 
and  to  study  surface  blast  phenomena.  The  variables  selected 
by  Schmidt  and  Holsapple  (1978)  were  apparent  crater  volume 
(V)  ,  the  heat  of  detonation  per  unit  mass  (Q) ,  initial  density 
of  the  explosive  (5) ,  initial  density  of  the  soil  (p)  , 
material  strength  of  the  soil  (Y),  spherical  explosive  charge 
radius  (a) ,  depth  of  charge  burial  (DOB) ,  and  the 

gravitational  constant  (G)  .  The  pi  terms  used  in  their 
analysis  are  listed  in  Table  4.1.  Crater  radius  (r)  and  depth 
(d)  are  also  analyzed  using  the  same  set  of  variables  such 
that : 


n,  =  n, 


fv 


n,  =  71, 


=  d(-P-) 


W 


(4.1) 


or 


TABLE  4.1.  CENTRIFUGE  CRATERING  PI  TERMS. 


I.  Surface  (Half-Buried)  Detonations 
(Schmidt  and  Holsapple,  1980) 


A. 

Mass  Set 

7ti  = 

I 

''  W 

(^)  (^) 

71-.  = 

C,  ft 

y 

7t„  = 

£ 

4 

ft 

Ciateiing  Efficiency 
Gravity  Scaled  Yield 


B .  Energy  Set 


71.  = 
nj  = 

71^  = 


C .  Gravity  Set 
71, 


71, 


71,  = 


E 

x 

£ 

ft 


V(^)  3/4 

E 

-L  I 

0.  ^ 

Y 


7t, 


£ 
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TABLE  4.1. 


CENTRIFUGE  CRATERING  PI  TERMS 
(CONTINUED) . 


II.  Buried  Detonations 

(Schmidt  and  Holsapple,  1977) 


Jii,  TCj,  from  I. A,  the  mass  set,  and 
^  W 


III.  Material  Strength  Model 

(Schmidt  and  Holsapple,  1979) 


Til,  71^  from  I. A,  the  mass  set,  and 


Tij  =  tan4> 


when  tests  of  half-buried  charges  are  performed  in  the 
same  soil  with  the  same  explosive,  713=0  and  both  and  7I5  are 
constant.  Thus,  using  the  Buckingham  Pi  theorem  (Buckingham, 
1914),  Schmidt  and  Holsapple  (1978)  wrote  the  functional 
relationship  as ; 

=  F  (4.2) 

Schmidt  and  Holsapple  (1978)  found  the  power  law  could  be 
used  to  relate  Ttj  to  7:,. 
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(4.3) 


=  k  -=  constant 

The  expansion  of  the  pi  terms  to  their  physical  definitions/ 
resulted  in  scaling  rules  for  dynamically  nonsimilar  tests  in 
Ottawa  sand; 

=  if  (i!^)  -»  {4.4a) 

0  Q  0 

r  =  if  (- (4.4b) 

^  P  0  6 

d  =  k.(-^)i/M- 

^9  0  o 

For  the  apparent  crater  dimensions,  Schmidt  and  Holsapple 
(1978)  defined  the  yield  exponents  as 


Pj^(1-0/3> 

(4.5a) 

-fi)  / 3 

(4.5b) 

/3 

(4 .50 

where  the  exponents  a,  P,  and  7  represent  the  material 
properties  of  the  soil. 

Using  the  power  law,  Schmidt  and  Holsapple  (1578)  found 
their  data  agreed  with  the  1  G  test  data  reported  by 
Piekutowski  (1974,  1975).  Figure  4.1  shows  the  agreement 
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LEGEND 


G/Q(W/i)l 


between  the  two  sets  of  data  for  apparent  crater  volume.  From 
such  figures,  Schmidt  and  Holsapple  (1978)  concluded  that 
large  explosive  yields  at  1  G  could  be  accurately  simulated  by 
small  charges  at  elevated  gravity,  and  that  for  half-buried 
charges  the  effect  of  increasing  the  gravitational 
acceleration  on  a  fixed  charge  is  equal  to  increasing  charge 
size  at  fixed  gravity. 

The  gravity-scaled  yield  (  (G/Q.)  (W/5)  accounted  for 
variation  in  explosive  properties  as  well  as  the  variation  in 
gravity.  Therefore,  Schmidt  and  Holsapple  (1978)  also 
concluded  tnat  use  of  a  gravity-scaled  yield  term  would  permit 
large  explosive  yields  to  be  simulated  by  a  small  scale 
equivalent  charge  of  a  different  type. 

When  all  other  variables  are  fixed,  the  power-dependence 
of  apparent  crater  dimensions  in  Otcawa  sand  on  charge  mass 
were  reported  by  Schmidt  and  Holsapple  (1978)  as  ^ for 
apparent  crater  volume,  ^  °  for  apparent  crater  radius, 
and  ®  for  apparent  crater  depth. 

Schmidt  and  Holsapple  (1990)  concluded  that  different 
scaling  laws  do  not  result  from  different  variables  and  pi 
term  sets.  The  scaling  laws  from  different  pi  term  sets 
appear  different  because  each  requires  different  testing 
procedures  to  hold  the  pi  terms  constant  and  to  maintain 
conditions  of  similarity.  If  all  of  the  variables  used  to 
relate  the  variables  in  one  pi  term  set  to  the  variables  in 
another  set  are  included  in  the  pi  term  sets,  any  one  set  of 


pi  terms  can  be  converted  to  any  other  through  the  product  of 
powers-transformation  procedure.  Thus,  as  long  as  specific 
energy  of  the  explosive  is  included  in  the  analysis,  charge 
mass  is  interchangeable  with  explosive  energy. 

If  the  soil  type  is  fixed,  a  variation  in  gravity  as 
required  by  similarity  produces  the  cube-root  scaling  law. 
The  quarter-root  scaling  law  will  only  result  if  gravity  is 
kept  constant,  but  no  two  similar  tests  {7I2,  Jtj,  =  constant) 
can  fulfill  these  restrictions.  Therefore,  Schmidt  and 
Holsapple  (1980)  showed  that  when  the  soil  type  is  fixed, 
quarter-root  scaling  is  impossible. 

2 .  MATERIAL  STRENGTH  MODEL 

Figure  4.2  (Holsapple  et  al.,  1978;  Schmidt  and 
Holsapple,  1979)  shows  that  when  the  cratering  efficiency  pi 
term  (nj  is  plotted  against  the  gravity  scaled  yield  pi  term 
(Ttj)  ,  those  materials  having  some  cohesion,  apparent  or  real, 
do  not  fall  on  the  curve  for  dry  sand.  At  low  G-leveJ.s  or  low 
explosive  yield,  the  cratering  efficiency  is  constant.  The 
apparent  crater  volume  is  proportional  to  the  charge  mass,  and 
the  cube  root  scaling  rule  applies. 

At  high  G-levels,  the  data  becomes  asymptotic  to  lines 
parallel  to  the  dry  sand  curve.  Beyond  some  G-level  or  charge 
size,  the  cube  root  scaling  rule  is  not  appropriate.  The 
crater  dimensions  exhibit  a  dependence  on  charge  mass. 
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LEGEND 


Various  Soil-Explosive  Combinati 
1978).  Tr„  =  Ve/W,  =  G/Q(W/A) 


To  account  for  this  observed  deviation  in  scaling, 
Schmidt  and  Holsapple  (1979)  used  an  energy  balance  approach 
to  derive  a  modified  gravity  scaled  yield  term: 

^2m  =  -•  (4.6) 

where  c  is  cohesion,  tan<))  is  the  frictional  shear  resistance, 
and  )c.i,  kj,  and  kj  are  dimensionless  constants.  first 
constant  (kj  was  defined  as  a  measure  of  th  j-k  done 
against  gravity.  The  second  constant  (kj)  was  defined  as  a 
measure  of  the  work  required  to  overcome  cohesion,  and  the 
third  constant  (k^)  was  defined  as  a  measure  of  the  work  done 
against  the  material  strength  increase  due  to  the  iiicrease  in 
confining  ptesauie.  Holsapple  et  al .  (1978)  found  ki 
empirically  for  Ottawa  sand  to  be  0.1.  Both  kj  and  kj  were 
found  to  be  approximately  1.0. 

Figure  4.3  is  a  graphical  presentation  of  the  material 
strength  model  for  crater  volume  (Holsapple  et  al .  ,  1978) .  At 
low  gravity  scaled  yields,  the  cohesion  dominates  the 
cratering  process.  The  curve  is  horizontal  and  crater  volume 
is  proportional  to  charge  mass.  The  cube  root  scaling  law 
applies.  Then  at  an  increased  gravity-scaled  yield  the  curve 
becomes  aligned  with  the  craters  formed  in  cohesionless  soil. 
The  soil  strength  due  to  the  frictional  shear  resistance 
dominates  the  cratering  process.  Holsapple  et  al .  (1978) 
concluded  their  analysis,  stating  that  the  observed  change  in 
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scaling  law  is  the  reason  small-scale  1  G  tests  in  soil  with 
cohesion  can  not  extrapolated  to  large-scale  1  G  tests. 

Holsapple  et  al .  (1978)  used  the  material  strength  model 
as  represented  by  to  explain  the  differences  in  crater 
volume  in  clay  and  in  sand.  Sand  has  a  relatively  high 
frictional  shear  resistance  and,  therefore,  has  a  strength 
that  increases  with  confining  pressure.  As  the  centrifugal 
acceleration  increases,  the  strength  of  the  sand  increases 
because  the  confining  pressure  increases.  In  the  short  term 
under  undrained  conditions,  clay  is  assumed  to  have  no 
frictional  shear  resistance,  but  a  relatively  large  cohesion. 
However,  cohesion  is  not  affected  by  gravitational 
acceleration,  so  the  strength  of  clay  remains  constant. 
Therefore,  Holsapple  et  al .  (1978)  stated  that  at  high  G- 
levels,  an  explosive  charge  will  produce  a  larger  crater  in 
clay  than  in  sand. 

Schmidt  and  Holsapple  (19/9)  stated  the  apparent  crater 
radius  and  depth  data  may  not  necessarily  fit  the  material 
strength  model  proposed  for  crater  volume.  Their  reason  was 
that  the  model  was  based  on  available  energy  and  the  work,  to 
excavate  the  soil.  No  parameters  for  shape  were  included. 
However,  the  radius  may  fit  the  model  to  some  extent  because 
volume  is  proportional  to  the  crater  depth  multiplied  by  the 
radius  squared. 

Overall,  Schmidt  and  Holsapple  (1979)  concluded  that 
subscale  craters  cannot  be  accurately  modeled  without 
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increased  gravity,  and  that  the  nsL-i.-. 
on  cohesion  and  frictional  shear  r»_oistance  accurately 
predicts  apparent  crater  volume  for  var'ous  soils  and 
explosive  sources. 

From  a  plot  of  scaled  crater  volume  (V/W)  versus  scaled- 
charge  mass  <N^W)  ,  Gaffney  et  al .  (1989)  concluded  that  the 

transition  point  between  strength  (cohesion)  dominated  scaling 
to  gravity  (frictional  shear  resistance)  dominated  scaling 
occurs  at  scaled  yields  of  a  few  kilograms.  However, 
differences  in  medium,  medium  density,  and  explosive  type  were 
ignored.  The  slope  of  the  decreasing  scaled  crater  volume  was 
found  to  be  -0.175  which  was  less  than  the  slope  of  -0.2 
calculated  by  Schmidt  et  al.  (1986). 

3 .  Subsurface  Detonations 

Schmidt  (i977)  used  four  pi  terms  to  analyze  charges 
buried  in  artist  modeling  clay  (Table  4.1).  Through  the 
application  of  the  product  of  powers  rule  which  permits  the 
combination  of  pi  terms,  a  scaling  law  for  crater  volume  in 
terms  of  the  depth  of  charge  burial  (DOB)  was  derived: 

^  Fh,nV^}  (4-7) 

The  scaling  law  was  shown  to  fit  the  data  from  86  cratering 
events  with  a  charge  yield  energy  range  of  12  orders  of 
magnitude.  Schmidt  and  Holsapple  (1977)  concluded  that  tests 
in  an  accelerated  reference  improve  the  prediction  of  large- 
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scale  events  and  attxibuted  the  accuracy  of  the  crater  volume 
scaling  law  to  the  gravity-scaled  yield  term  (tTj)  . 

Schmidt  (1979)  repeated  the  analysis  of  the  effects  of 
charge  burial  on  apparent  crater  volume.  When  the  gravity 
scaled  yield  (Ttj)  was  expanded  to  include  the  material 
strength  effects  of  cohesion  and  frictional  shear  resistance, 
the  scaling  law  became: 

=  Fin,)  (4.8) 

The  exponential  coefficient  for  on  the  right-hand  side  of 
Equation  (4.7)  did  not  vary  within  the  range  of  the  data,  and 
was  therefore  set  to  zero. 

B.  SLOPE  STABILITY  IN  CRATERING 

Schmidt  and  Holsapple  (1979)  added  that  crater  shape  will 
be  effected  by  slope  stability.  When  the  side  of  the  crater 
becomes  unstable,  slumping  will  occur.  Slope  suability  and 
slumping  are  functions  of  the  crater  depth,  cohesion,  and 
angle  of  internal  friction.  Therefore,  Schmidt  and  Holsapple 
(1979)  concluded  limits  exist  on  apparent  crater  radius  and 
depth  combinations.  The  lower  limit  on  the  aspect  ratio 
( radius/depth)  of  the  apparent  crater  was  defined  as  l/tan0 
for  desert  alluvium  where  <{)  is  the  angle  between  the  crater 
wall  and  the  original  ground  surface.  Cohesion  decreases 
crater  aspect  ratio  for  a  given  crater  depth. 
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c. 


MOISTURE  CONTENT  EFFECT  ON  CRATER  DIMENSIONS 


Schmidt  and  Holsapple  (1979)  also  analyzed  the  effect  of 
moisture  content  on  apparent  crater  dimensions  in  Desert 
Alluvium.  They  found  small  differences  in  moisture  content 
produced  large  differences  in  crater  size  and  structure.  A 
decrease  in  moisture  content  resulted  in  flatter  craters, 
whereas  at  a  larger  moisture  content  the  crater  snape  became 
more  hemispherical.  Tnus,  Schmidt  and  Holsapple  (1979) 
concluded  the  decrease  in  cohesion  at  the  lower  moisture 
content  caused  the  flattening  of  the  apparent  craters. 

D.  SCALE  EFFECT  IN  CRAfERING 

Goodings  et  al.  (1988)  concluded  there  are  no  detectable 
scale  effects  due  to  centripetal  acceleration  between  31  G  and 
101  G,  Given  the  same  soil  and  explosive,  craters  formed  on 
a  centrifuge  at  G-levels  ranging  from  1  G  to  463  G  vrill  all 
follow  the  same  scaling  law.  The  power  dependence  (yield 
exponent)  of  crater  dimensions  on  a  half-buried  charge  mass  in 
Ottawa  sand  were  found  to  be  w®®’  for  apparent  crater  volume, 
W'  ’’  for  apparent  crater  radius,  and  W“  for  apparent  crater 
depth . 

E .  SOIL  DENSITY  EFFECT  ON  CRATER  VOLUME 

Goodings  et  al .  (1980)  included  tests  at  different  soil 
densities  ranging  from  1590  Itg/m’  to  1720  kg/in^.  The  tests 
showed  an  8  percent  decrease  in  density  (36  pe’"C9nt  decrease 
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in  relative  density)  produced  a  72  percent  increase  in 
apparent  crater  volume. 


F .  CORIOLIS  EFFECT  ON  CRATER  SHAPE 

Kutter  et  al .  (1985)  studied  the  effect  of  blast  loading 
on  tunnels.  When  compass  directions  are  assigned  to  the 
specimen  surface  (E'igure  4.4),  the  authors  found  the  apparent 
crater  profiles  exhibited  "marked  asymmetry"  in  the  East-West 
direction,  and  attributed  this  distortion  to  the  Coriolis 
effect  and  the  close  proximity  of  the  tunnel.  The  apparent 
crater  profile  in  the  North-South  direction  remained 
symmetrical . 

Schmidt  et  al .  (1986)  and  Steedman  (1990)  analyzed  the 
Coriolis  effect  on  the  flight  of  the  centrifuge  crater  ejecta. 
In  both  analyzes,  particle  flight  was  a  function  of  ejection 
angle  and  particle  velocity.  Schmidt  et  ai  .  (1986)  stated  the 
Coriolis  acceleration  shortens  the  ballistic  range  of  the 
particles  ejected  towards  the  East  and  extends  the  ballistic 
range  of  the  particles  ejected  towards  the  West.  Thus  within 
the  crater,  the  ejecta  will  become  more  concentrated  in  the 
West.  The  result  is  an  asymmetrical  apparent  crater  profile. 

Schmidt  et  ai .  (1986)  observed  asymmetrical  crater 
profiles  in  dry  Ottawa  Flintshot  Sand  where  the  ejection 
angles  typically  ranged  from  40  to  45  degrees.  The  Coriolis 
effect  was  shown  to  increase  with  ejection  angle.  However, 
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Centrifuge 

Rotation 


Figure  4 


.4  Compass  Directions  Assigned  to 

the  Centrifuge  Specimen  as  Pictured 
from  the  Center  of  Rotation. 
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Schmidt  et  al .  (1986)  stated  that  below  55  to  60  degrees,  the 
Coriolis  effect  should  not  be  significant. 

Steedman  (1990)  assumed  the  trajectory  of  the  ejecta  was 
a  straight  line  directed  ahead  of  the  centrifuge  bucket  due  to 
the  initial  tangential  velocity  (Figure  4.5).  Thus,  at 
certain  combinations  of  centrifuge  speed  and  particle 
velocity,  the  particles  will  return  to  the  bucket.  Steedman 
(1990)  attributed  the  accumulation  of  ejecta  ahead  of  the 
original  detonator  location  to  the  high  absolute  velocity  of 
the  particles  and  the  shorter  distance  traveled  by  the 
particles  relative  to  the  centrifuge  bucket. 

Steedman  (1990)  noted  that  his  analysis  did  not  include 
the  strong  wind  created  by  the  rotation  of  the  centrifuge 
which  may  prevent  the  ejecta  from  following  the  theoretical 
straight  line  trajectories.  From  a  series  of  centrifuge  tests 
designed  to  model  the  effects  of  blast  loading  on  piles 
located  in  saturated  sand,  Steedman  concluded  there  was  no 
evidence  that  the  Coriolis  effect  distorted  the  shape  of  the 
apparent  crater. 

Figure  4.6  shows  the  apparent  crater  profiles  measured  by 
Gill  and  Kuennen  (1991)  in  the  East-West  direction.  The 
average  crater  profile  (test  number  500-00  ave . )  of  the  tests 
produced  by  a  half-buried  charge  is  symmetrical  about  a 
vertical  axis  located  at  the  origin  of  the  blast.  In 
contrast,  the  crater  profiles  produced  by  buried  charges  all 
appear  to  shift  to  the  right.  The  centiifuge  was  reported  by 
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Figure  4.6  Scaled  Crater  Profiles  Measured  by  Gill 
and  Kuennen  (1991). 
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Gill  and  Kuennen  to  rotate  clockwise  so  that  the  movement  of 
the  specimen  was  from  right  to  left.  Thus,  the  apparent 
crater  profile  from  buried  charges  appear  to  become  skewed  in 
the  direction  opposite  the  direction  of  centrifuge  rotation. 
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SECTION  V 


MODEL  DESCRIPTION 

A.  EXPLOSIVES 

Explosive3  used  in  this  research  were  RP-83  detonators 
specially  prepared  by  Reynolds  Industries  Systems,  Inc  of  San 
Ramon,  California.  Four  different  sizes  of  the  detonator  were 
used  to  model  several  explosive  weights.  The  make-up  of  each 
detonator  is  shown  in  Table  5.1.  The  output  charge  consists 
of  a  multiple  of  22C  mg  cyclotetramethl  enetetranitrarriine  (PBX 
9407)  pres.sings  (Figure  5.1).  Each  220  mg  pressing  is  4.4  mm 
in  length  by  7.0  mm  in  diameter.  Thus,  each  larger  detonator 
is  4.4  mm  longer  than  the  next. 

The  total  charge  mass  is  converted  to  an  equivalent 
weight  of  TNT  for  the  data  analysis.  PBX  9407  is  94  percent 
cyclotrimechylene  trinitramine  (RDX)  and  6  percent  Exon  461 
(Baker  et  al . ,  1980)  .  Based  on  the  heat  of  detonation,  the 
TNT  equivalent  weights,  for  PBX  9407  and  RDX  are  1.136  and 
1.149  respectively.  For  cratering,  the  Air  Force  uses  a  factor 
of  1.115  to  convert  RDX  to  TNT.  Thus,  the  equivalent  mass  of 
TNT  can  be  expressed  as: 


(1.09)  PBX 


TNT . 


(5.1) 
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B. 


TEST  MATRICES 


Table  5.2,  shows  the  test  matrix  for  this  research.  The 
G-levels  (1  G  =  9.81  m/s^)  presented  in  the  matrix  are  the  G- 
levels  at  the  center  of  mass  of  the  detonator. 

1 .  Constant  Normalized  Depth  of  Burial. 

a  ^  0.73  m/ka‘^^) 

Detonators  were  buried  at  a  normalized  depth  of 
charge  burial  (X)  equal  to  0.73  m/kg^^^  to  insure  as  much  of 
the  explosive  energy  as  possible  was  released  into  the  soil. 
Drake  and  Little  (1983)  introduced  an  equivalent  effect 
coupling  factor  to  adjust  the  explosive  energy  directed  into 
a  medium  according  to  the  depth  of  charge  burial.  For  soil, 
the  maximum  explosive  energy  is  shown  to  be  released  into  the 
medium  at  normalized  depths  of  burial  of  0.56  m/kg^'^  and 
greater. 

The  centrifuge  test  series  was  designed  to  study 
explosive-induced  stress  wave  transmission  and  crater  size  as 
a  function  of  saturation.  This  test  series  formed  the  first 
part  of  the  matrix  (Table  5.2a)  .  The  modeling  of  models 
technique  was  used  to  simulate  a  prototype  charge  mass  of  7  kg 
of  TNT.  The  largest  (880  mg)  RP-83  detonators  were  fired  at 
19  G  and  the  smallest  (220  mg)  RP-83  detonators  were  fired  at 
26  G.  Each  detonator  was  fired  in  dry  specimens  and  spec.i.mens 
saturated  at  compaction  levels  of  20,  40,  60,  and  70  percent. 


TABLE  5.2a.  TEST  MATRIX  FOR  THE  CONSTANT 

NORMALIZED  DOB  SERIES  (k  =  0.73  m/kg^'^) 


Two  more  test  series  were  added  to  include  tests  at 
higher  G-levels.  The  series  modeled  49  and  124  kg  of  TNT,  and 
subjected  the  specimens  to  49  G  and  67  G  respectively.  In 
these  series,  the  tests  were  conducted  at  compaction 
saturation  levels  of  0,  20,  40,  and  60  percent 

Last,  six  tests  were  conducted  without  centrifugation  at 
1  G  (Table  5,2b)  .  Tests  were  conducted  at  1  G  because  1  G 
tests  aie  legitimate  tests  in  and  of  themselves,  and  to  show 
how  the  use  of  a  centrifuge  changes  the  results.  From  the 
centrifuge  scaling  relations  it  is  apparent  the  1  G  tests  are 
in  no  way  related  to  the  centrifuge  tests.  The  1  G  tests  only 
model  a  detonator  buried  in  a  few  centimeters  of  sand,  whereas 
the  centrifuge  tests  model  several  kilograms  of  explosive 
buried  in  a  few  meters  of  sand. 

The  compaction  saturation  levels  selected  were  0,  20,  and 
60  percent.  All  of  the  tests  were  conducted  in  the  centrifuge 
buckets.  Two  of  the  1  G  tests,  one  at  20  percent  compaction 
saturation  and  one  at  60  percent  compaction  saturation  were 
spun  on  the  centrifuge  to  57.3  G  before  being  tested.  The  two 
1  G  tests  conducted  in  dry  beach  sand  did  not  have  the  same 
method  of  compaction.  One  specimen  was  vibrated  while  the 
other  was  pluviated. 


2 .  Varied  Normalized  Depth  of  Burial 
(X  =  0.23  to  0.99 

A  separate  series  of  11  tests  was  conducted  to  study 
the  variation  in  crater  size  and  shape  with  depth  of  charge 
bur.1  1  (Table  5.2c).  In  this  test  series,  X  ranged  from  0.23 
m/kg*^^  to  0.99  m/kg^^^.  All  of  the  tests  used  oven-dry  Tyndall 
Beach  Sand,  and  all  of  the  specimens  were  compacted  through 
plaviation.  In  the  design  of  this  test  matrix,  the  modeling 
of  models  techiiique  was  used  whenever  possible. 


SECTION  VI 


SPECIMEN  PREPARATION 
A.  SOIL  DESCRIPTION 

The  soil  used  in  this  research,  termed  Tyndall  Beach 
Sand,  was  quartz  beach  sand  found  on  Tyndall  Air  Force  Base  in 
Florida.  The  sand  was  washed,  oven-dried,  and  passed  through 
a  number  20  sieve  prior  to  use.  As  shown  in  Figure  6.1,  this 
sand  has  an  almost  uniform  particle  size  distribution.  The 
characteristic  shape  of  the  sand  is  classified  as  subrounded 
to  subangular  (Figure  6.2).  The  calculated  parameters  for 
Tyndall  Beach  Sand  are  listed  in  Table  6.1. 


TABLE  6.1.  SOIL  PARAMETERS  FOR  TYNDALL  BEACH  SAND. 


The  sand  was  compacted  either  through  vibration  or 
pluviation  to  a  dry  density  of  1522  kg/m^.  The  minimum  and 


mdall  Bea 


Figure  6.2  Photograph  Showing  the 

Characteristic  Shape  of  Tyndall 
Beach  Sand. 
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maximum  dry  densities  for  the  sand  are  1450  kg/m’  and  1630 
kg/m^  respectively  (ASTM  D4253  and  D4254) .  At  a  dry  density 
of  1522  kg/m^,  the  void  ratio  is  0.74  which  leads  to  a 
relative  density  of  39  percent. 

B.  SATURATION  PROCEDURE 

For  tests  in  partially  saturated  sand,  distilled  water 
was  added  to  the  specimen  prior  to  placement.  The  procedure 
involved  weighing  out  enough  dry  sand  for  one  lift  in  the 
centrifuge  buckets,  adding  the  desired  weight  of  distilled 
water,  and  mixing  the  moist  sand  in  a  Chiyoda  flexible  wall 
mixer  (Omni-Mixer  model  OM-10,  Figure  6.3). 

Five  lifts  were  used  to  fill  the  aluminum  centrifuge 
buckets  measuring  457  mm  in  inside  diameter  and  305  mm  in 
height.  The  first  four  lifts  were  50  mm  in  height.  The 
height  of  the  fifth  lift  was  determined  by  the  normalized 
depth  of  charge  burial  (X.)  which  in  turn  was  determined  by  the 
test  G-level  and  charge  size.  Therefore,  for  similar  tests, 
the  depth  of  the  overburden  was  varied  instead  of  the  charge 
elevation . 

C  .  PLACEMENT  METHODS 

All  partially  saturated  specimens  were  compacted  using 
the  vibratory  method.  In  the  constant  normalized  DOB  test 
series,  the  dry  specimens  were  compacted  using  the  vibratory 
method  to  maintain  a  consistent  method  of  preparation.  In  the 
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test  series  spanning  a  range  of  DOB,  all  specimens  were 
pluvitated. 

1 .  Vibration 

A  top-down,  vertical  vibration  technique  was  used  to 
compact  the  partially  saturated  specimens.  A  Syntron  magnetic 
vibrator  (model  No.  V51  Dl)  was  attached  to  an  circular 
aluminum  plate  cut  to  0.44  m  in  diameter.  Each  lift  was 
compacted  separately  by  placing  the  plate  down  against  the 
sand  (Figure  6.4). 

The  compaction  was  controlled  by  lift  height 
markings  painted  around  the  inside  of  the  bucket.  Once  the 
leveled  lift  matched  the  appropriate  markings,  the  dry  density 
of  1522  kg/m^  had  been  attained.  As  a  final  check,  the  full 
bucket  was  weighed  at  the  completion  of  compaction. 

2 .  Pluviation 

The  pluviator  used  in  this  research  was  designed  by 
Dr.  Taylor  of  Virginia  Polytechnic  Institute  and  State 
University^  and  is  shown  in  Figure  6.5.  In  pluviation,  the 
soil  is  rained,  allowed  to  free-fall,  through  holes  and 
several  sieves  (wire  screens)  into  the  specimen  bucket  (Figure 
6.6).  For  a  specific  soil  and  set  free-fall  heights,  the 
firial  density  of  the  soil  as  it  falls  into  the  bucket  is  a 

^Personal  Communication,  Mike  Purcell,  Engineering  Assistant, 
Applied  Research  As  i  af  es ,  ,  'T’yndall  AFB,  Florida. 
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Figure  6.4  FMC  Corporation  Syntron  Magnetic 
Vibrator . 
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function  of  the  number  and  size  of  the  holes  in  the  board,  and 
the  opening  sizes  of  the  sieves  (Fragaszy  and  Taylor,  1989) . 
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SECTION  VII 


TEST  PROCEDURES 
A.  CENTRIFUGE  FACILITY 

The  centrifuge  facility  at  Tyndall  Air  Force  Base  was 
utilized  for  this  research.  The  facility  contains  a  133  G-KN 
centrifuge  (Genisco  model  E-185)  enclosed  in  a  circular 
concrete  structure.  The  centrifuge  radius  is  1.83  m,  and  the 
payload  platforms  each  measure  7  6  cm  x  7  6  cm.  The  payload 
platforms  are  attached  to  the  centrifuge  by  two  symmetrical 
cantilever  arms  and  are  thus  allowed  to  swing  freely  during 
the  operation  of  the  centrifuge  (Figure  3.1.b).  Weights  were 
added  to  the  opposing  arm  to  balance  the  centrifuge.  An 
electronic  balancing  system  assisted  this  manual  balancing 
process  and  automatically  kept  the  centrifuge  in  balance 
during  flight. 

The  data  acquisition  system  (Pacific  Transient  Data 
Recorder  model  5700)  and  a  Miletus  Associates  (model  0)  video 
camera  have  been  mounted  on  the  central  axis  of  the  centrifuge 
(Figure  7.1).  These  additions  allow  on-board  data  capture  and 
the  photographic  recording  of  all  centrifuge  tests.  The 
centrifuge  is  operated  by  a  hydraulic  drive  system  from  a 
nearby  control  room.  The  video  is  connected  to  a  monitor  in 


Figure  7.1  Miletus  Video  Camera  and  Pacific 

Data  Acquisition  System  Mounted  on 
the  Centrifuge. 
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the  control  room  such  that  real-time  viewing  of  a  test  is 
possible . 

B,  hetonator  placement 

The  RP-83  detonators  were  placed  in  the  specimens 
just  prior  to  centrifuge  start-up.  Detonator  placement  could 
not  be  made  during  compaction  because  of  safety  regulations. 

For  the  partis  saturated  specimens,  an  8  mm  diameter 
hole  was  manually  drilled  into  the  sand.  Capillary  forces 
within  the  sand  were  sttoi;g  enough  to  prevent  the  walls  of  the 
hole  from  caving  in.  For  the  dry  specimens,  a  glass  tube  with 
an  inside  diameter  of  10  mm  was  used  to  retain  the  surround4.:ig 
soil  and  permit  excavation.  U.sing  a  guide  and  suction  to 
remove  sand  grains  from  inside  the  glass  tube,  the  glass  tube 
was  gently  pushed  into  t)ie  sand. 

When  the  detonator  was  placed  into  the  drilled  hole  or 
the  glass  tube,  a  mark  on  the  detonator  wires  corresponding  to 
the  sand  surface  verified  the  correct  placement  of  the 
detonator.  The  guide  designed  to  place  the  detonator  at  the 
centt...  or  the  bucket  and  to  keep  the  drill  or  glass  tube 
oriented  in  the  vortical  direction  is  pictured  in  Figure  7.2. 
Markings  along  the  side  of  the  guide  a] lowed  the  depth  to  the 
botton  of  the  hole  to  be  measured  within  1  mm.  All  depth  of 
burial  dis'aiices  are  measured  Irorn  ttie  surxace  to  the  center 
of  muss  of  the  charge. 


0  7 


70 


Dry  sand  was  poured  back  into  the  hole  to  keep  the 
detonator  vertical  and  to  hold  the  detonator  in  place. 
Knudson  et  al .  (1972)  found  that  the  choice  of  material  used 
to  stem,  backfill,  a  charge  was  not  critical,  and 
that significant  increases  in  crater  radius  or  depth  did  not 
occur  with  stemiaing  beyond  50  percent  of  the  emplacement  hole 
length.  For  the  dry  tests,  the  glass  tube  was  removed  once 
the  hole  had  been  stemmed. 

C .  CENTRIFUGE  OPERATION 

The  operation  of  the  cencrifuge  du’^ing  a  test  consisted 
of  running  the  ce.ntrifuge  up  to  the  desired  G-level,  firing 
the  detonator,  and  decelerating  the  centrifuge  back  dov/n  to  a 
statio.nary  position.  Each  test  required  approximately  10 
minutes  from  centrifuge  start-up  to  shut-down. 

The  electronic  connection  of  the  Reynolds  Industries  EBW 
firing  system  (model  FS-17)  to  the  detonator  passed  through 
one  of  the  28  slip  rings  of  the  centrifuge.  The  firing 
sequence  was  manually  armed  and  operated.  When  the  desired  G- 
level  was  attained,  the  charge  was  detonated. 

D  .  APPAFENT  CRATER  MEA.nURF.MENT 

1 .  Inside  Diameter  and  Depth 

The  crater  depth  and  inside  diameter  were  measured 
immediately  after  the  centrifuge  came  to  rest.  No  changes  in 
the  apparent  crater  shape  were  observed  daring  centrifuge 


deceleration.  The  measurements  were  taken  in  two  different 
directions,  one  perpendicular  to  the  other.  Thus,  the 
reported  dimensions  are  an  average.  Both  depth  measurements 
were  taken  at  the  deepest  point  in  the  crater. 

2 .  Volume 

Crater  volume  measurements  were  made  from  a  clay 
replica.  The  clay  replica  was  made  from  a  concrete  mold  taken 
after  the  crater  diameter  and  depth  had  oeen  measured.  When 
the  clay  replica  was  filled  with  water,  the  weight  of  water  in 
grams  equaled  the  crater  volume  in  cubic  centimeters, 

E.  HOISTURE  CONTENT  MEASUREMENT 

The  downward  migration  of  water  that  occurs  during  a 
centrifuge  test  quickly  became  evident.  The  centripetal 
acceleration  of  a  centrifuge  desaturates  the  specimeii  as 
gravity  would  desaturate  the  prototype  soil  layer.  To 
determine  the  extent  of  the  desaturation,  the  specimen  was 
sampled  inamediately  following  the  test .  The  moisture  content 
present  at  the  completion  of  the  test  was  assumed  ro 
approximate  the  moisture  content  present  during  the  test,  and 
the  upward  migration  of  water  during  the  deceleration  of  the 
centrifuge  was  assumed  to  be  minimal. 

Samples  were  taken  from  the  top,  middle,  and  bottom  51  rrmn 
of  the  specimen.  After  being  weighed,  the  samples  were  dried 
in  an  oven  at  120  °C  for  24  hours.  One  hole  was  located 
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towards  the  center  of  the  bucket  at  the  crater  edge,  while  a 
second  was  located  at  the  side  of  the  bucket.  Initially,  four 
holes  were  used  to  sample  the  specimen,  two  at  the  12  o'clock 
position,  and  two  at  the  3  o'clock  position  (Figure  7,3)  . 
After  several  tests,  however,  four  holes  proved  to  be 
redundant.  The  moisture  content  at  the  crater  edge  or  at  the 
side  of  the  bucket  did  not  vary  significantly  from  one 
position  to  the  other;  therefore,  sampling  at  the  3  o'  clock 
position  was  discontinued. 
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SECTION  VIII 


RESULTS  AND  ANALYSIS 

A.  SCALED  CRATER  DIMENSION  VERSUS  SCALED  EXPLOSIVE  MASS 

The  current  research  differs  from  the  previous  centrifuge 
cratering  studies  in  soil  type,  saturation,  and  density  (Table 
8.1).  Figures  8.1,  8.2,  and  8.3  show  the  scaled  apparent 
crater  volume,  radius,  and  depth  plotted  against  scaled 
explosive  mass.^  The  best-fit  lines  for  these  relationships 
developed  from  the  current  research  have  the  same  general 
slope  and  plot  above  those  developed  from  the  data  presented 
by  Goodings  et  al.  (1988)  and  Schmidt  and  Holsapple  (1978, 
1980)  .  The  constant  slope  confirms  the  independence  of  soil 
type,  saturation,  and  density  from  explosive  mass. 

Schmidt  and  Holsapple  (1980)  and  Goodings  et  al .  (1988) 
were  able  to  relate  1  G  tests  with  ceiitri.fuge  tests  because 
the  charges  were  placed  at  the  soil  surface  (half-buried) . 
Surface  charges  have  confining  pressures  of  zero  at  all  G- 
levels.  For  buried  charges,  the  confining  pressures  are  a 
function  of  depth  and  G-level  which  makes  the  comparison  of 
tests  at  1  G  and  high-G  Uiireaoonable  for  frictional  materials 


■’Plotted  in  english  and  metric  units  to  allow  comparison  with 
other  work,  conversion  equations  provided  in  figures.  Scaled 
crater  dirr.en.sions  are  prototype  crater  dimensions. 


TABLE  8.1.  CENTRIFUGE  TEST  SPECIFICATIONS. 


Current 

Goodings 

Schmidt  & 

Schmidt  & 

Research 

et  al . 

Holsapple 

Holsapple 

(1992) 

(1988) 

(1978) 

(1979) 

Soil 

Tyndall 

Ottawa 

Ottawa 

Desert 

Type 

Beach 

Sand 

Sand 

Sand 

Alluvium 

Sat. (%) 

0-70 

0 

0 

15  " 

Dry 

Density 

1522 

1704 

1779 

1572  = 

(kg/m3) 

DOB 

Fully 

Half- 

Half- 

Half- 

(m/kg^'^) 

Buried 

Buried 

Buried 

Buried 

(0 . 23 

(0) 

(0) 

(0) 

to  0.99) 

a.  Average  value  for  21  tests  with  saturations  ranging 
from  3  to  29  percent. 

b.  Average  value  for  19  tests  with  dry  densities  ranging 
from  1677  kg/m’  to  1802  kg/m^. 

c.  Average  value  for  21  tests  with  dry  densities  ranging 
from  1429  kg/m^  to  1657  kg/m\ 
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aled  Apparent  Crater  Volume  Versus  Scaled  Explosive  Mass 


Scaled  Explosive  Mass,  TNT  Equivalent  (kg) 
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Figure  8.2  Scaled  Apparent  Crater  Radius  Versus  Scaled  Explosive  Mass, 

TNT  Equivalent.  (Conversions:  0.305ft  =  1.0m,  0.4541b  =  1.0kg) 


Explosive  Moss,  TNT  Equivolent  (kg) 


(lu)  M;d9Q  ja;DJ3  ^uajoddy  p9|Dos 


Scaled  Apparent  Crater  Depth  Versus  Scaled  Explosive  Mass# 

TNT  Equivalent.  (Conversions:  0.305ft  =  1.0m,  0.4541b  =  l.Okg) 


including  dry  sand.  For  cohesive  materials,  a  comparison  of 
1  G  and  high-G  tests  may  not  be  meaningless.  As  shown  in 
section  8.B,  the  capillary  forces  within  the  sand  produce  an 
apparent  cohesion  which  controls  crater  formation  at  elevated 
G-levels  as  well  as  at  1  G.  Therefore,  the  partially 
saturated  1  G  tests  were  included  in  the  derivation  of  the 
centrifuge  scaling  laws  whereas  the  dry  1  G  tests  were  not. 

The  partially  saturated  tests  are  also  compared  to  tests 
in  moist  (=  4  percent  moisture  content,  =  15  percent 

saturation)  desert  alluvium  (Schmidt  and  Holsapple,  1979) 
(Figures  8.4,  8.5,  8.6).  Again,  the  best-fit  lines  for  the 

data  from  the  current  research  are  greater  than  and  parallel 
to  the  best-fit  lines  of  other  centrifuge  research.  The  best- 
fit  lines  drawn  through  the  data  on  a  log-log  scale  have  the 
fom  Y  O'  where  the  slope  n,  in  these  figures,  is  the  yield 
exponent  (see  section  II),  VJ  is  the  explosive  mass,  and  Y  is 
the  scaled  apparent  crater  dimension  (volume,  radius,  depth)  . 

Table  8.2  lists  the  yield  exponents  from  this  research  as 
well  as  from  the  research  by  Goodings  et  al .  (1988),  Schmidt 
and  Holsapple  (1978,  1980),  and  the  U.S.  Army  Engineer 

Waterways  Experiment  Station  (Strange  et  al .  ,  1961  ;  Rooke  et 
al.,  1974) .  The  yield  exponents  of  the  current  research  are 
closely  aligned  with  the  yield  exponents  of  the  other 
research.  The  yield  exponents  for  the  dry  Tyndall  Beach  Sand 
are  almost  identical  to  the  yield  exponents  of  the  dry  Ottawa 
Sand,  and  the  yield  exponents  for  the  moist  Tyndall  Beach  Sand 
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Scaled  Explosive  Mass,  TNT  Equivalent  (kg) 


Scaled  Apparent  Crater  Radius  Versus  Scaled  Explosive  Mass, 

TNT  Equivalent.  (Conversions:  0.305ft  =  1.0m,  0.4541b  =  1.0kg) 
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are  in  agreement  with  the  yield  exponents  for  the  moist  Ottawa 
Sand  and  moist  Desert  Alluvium.  In  contrast  to  the  WES  report 
(Strange,  1961),  apparent  crater  radius  and  depth  both  appear 
to  scale  by  0.3  rather  than  by  1/3. 

The  specimens  tested  in  the  current  research  were 
compacted  to  a  dry  density  of  1522  kg/m^  which  is  less  than 
the  1704  kg/m^  used  by  Goodings  et  al .  (1988)  and  less  than 
the  1779  kg/m^  for  Ottawa  Sand  used  by  Schmidt  and  Holsapple 
(1978,  1980)  (Table  8.1).  Therefore,  the  data  plot  as 

expected,  above  the  data  reported  by  Goodings  et  al.  (1988) 
and  Schmidt  and  Holsapple  (1979,  1980) . 

Table  8.3  lists  by  depth  of  charge  burial  (DOB),  the 
yield  exponents  for  dry  to  moist  sand,  "wet"  (unsaturated) 
sand,  and  moist  desert  alluvium  established  by  Rooke  et  al. 
(1974)  for  full-scale  tests.  No  consistent  trend  can  be 
identified  in  these  data  for  a  variation  in  yield  exponent 
with  DOB.  However,  the  yield  exponents  generated  by  the 
current  research  are  in  reasonable  agreement  with  the  variable 
results  reported  by  Rooke  et  al.  (1974). 

The  correlation  coefficients  (R^  values)  from  a 
regression  analysis  using  the  method  of  least  squares  on  the 
log  of  the  current  research  data  for  the  constant  DOB  series 
(X  =  0.73  m/kg^'^)  are  listed  in  Table  8.4.  The  R^  values  for 
the  crater  dimensions  in  dry  sand  are  lower  than  the  values 
for  the  crater  dimensions  in  partially  saturated  sand  in  part 
because  no  1  G  tests  are  included.  All  of  the  R^  values  for 
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TABLE  8.3.  YIELD  EXPONENTS  REPORTED  BY  WES  (Rooke  et 
al.  ,  1974)  . 


Normalized  DOB 

Crater 

Dimension 

Dry  to 
Moist 

Wet 
Unsat 
Sand  Sand 

Moist 

Desert 

Alluvium 

0.05  <  X  <  0.2 

Radius 

Depth 

Volume 

0.280 

0.270 

0.906 

0.304 

0.313 

1 . 012 

0.2  <  X  <  0.5 

Radius 

Depth 

Volume 

0.279 

0.262 

0.815 

0.323 

0.263 

0.5  <  X  <  0.9 

Radius 

Depth 

Volume 

0.365 

0.266 

1 . 059 

0 . 259 
0.327 

0 . 840 

0.9  <  X  <  1.1 

Radius 

Depth 

Volume 

0.268 

0.237 

0 .282 
0.280 
0.859 

1.1  <  X  <  2.0 

Radius 

Depth 

Volume 

0.305 

0.238 

0.863 

0.306 

0.335 

0.951 

Table  8.4.  CORRELATION  COEFFICIENTS  FOR  THE  CURRENT 
RESEARCH  AT  ALL  G-LEVELS . 


Dry 

Partially 

Saturated 

Volume 

=  82.5% 

Volume 

=  99.8% 

Radius 

R^  =  93.9% 

Radius 

r2  =  99.8% 

Depth 

R^  =  62.3% 

Depth 

=  98.2% 
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the  crater  dimensions  in  partially  saturated  sand  are  above 
0.95.  The  lower  observed  scatter  in  the  data  presented  by 
Goodings  et  al .  (1988)  may  be  attributed  to  a  higher  soil 
density.  Greater  consistency  is  expected  for  denser 
specimens . 

B.  MATERIAL  STRENGTH  MODEL 

The  current  research  data  at  compaction  saturation  (S^) 
levels  of  20  and  60  percent  are  applied  to  the  material 
strength  model  proposed  by  Schmidt  and  Holsapple  (1979)  (see 
section  III).  Figures  8.7,  and  8.8  show  the  pi  terms  for 
apparent  cratering  efficiency  (K,  =  Vp/W) ,  apparent  crater 
radius  (tc^  =  r(p/W)^''^)  and  apparent  crater  depth  (71^  = 
d(p/W)^'^)  plotted  against  the  gravity  scaled  yield  pi  tern  (Ttj 
=  G/Q(W/5]  ^^^)  .  These  figures  follow  the  material  strength 
model  presented  in  Figure  4.3.  In  the  derivation  of  the 
gravity  scaled  yield,  the  soil  density,  explosive  density,  and 
explosive  specific  energy  were  ta)cen  to  be  1522  )cg/m\  1467 
hg/m^,  and  5.35  x  10^°  ergs/g  (Meyer,  1987)  respectively. 

At  60  percent  compaction  saturation,  the  cratering 
efficiency  appears  to  follow  the  model  closely  (Figure  8.7.b) . 
At  low  yields,  the  curve  is  horizontal  and  the  apparent 
cohesion  dominates  the  cratering  process.  Then  at  a  gravity 
scaled  yield  of  approximately  3  x  10"\  the  frictional 
component  of  soil  strength  begins  to  dominate  the  cratering 
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Figure  8,7  Cratering  Efficiency  (fTv )  Versus  Gravity 
Ccaled  Yield  (1^2)/  Tyndall  Beach  Sand: 
a)  S  =  0%,  20%;  b)  S  =  0%,  60%. 


88 


Crater  Radius,  rfi  =  Tfr  Crater  Depth,  lTi 
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process.  The  apparent  crater  volume  declines  with  increasing 
gravity  scaled  yield. 

Figure  8. 7. a  showing  cratering  efficiency  at  20  percent 
saturation  is  not  as  definitive.  The  transition  point  between 
cohesion  and  friction-dominated  cratering  appears  to  be  just 
at  or  beyond  the  yields  tested.  Thus,  there  is  an  indication 
that  cohesion  dominates  the  cratering  process  to  a  higher 
yield  at  20  percent  compaction  saturation  than  at  60  percent 
compaction  saturation.  The  additional  moisture  in  the 
specimens  compacted  at  60  percent  may  have  a  weakening  effect. 
Overall,  the  cratering  efficiency  data  support  the  results  of 
Gaffney  et  al .  (1989).  In  Tyndall  Beach  Sand,  the  transition 
point  occurs  at  a  few  kilograms  of  explosive  not  a  few 
thousand  kilograms  as  reported  by  Schmidt  et  al .  (1986). 

Apparent  crater  depth  (Figure  8. 8. a)  also  appears  to  follow 
the  material  strength  model  proposed  by  Schmidt  and  Holsapple 
(1979)  which,  for  reasons  mentioned  in  section  IV,  is 
unexpected.  As  for  the  cratering  efficiency,  the  transition 
yield  point  for  apparent  crater  depth  at  60  percent  compaction 
saturation  is  below  the  transition  yield  point  at  20  percent 
compaction  saturation. 

Figure  3.8.b  showing  apparent  crater  radius  is  less  clear 
because  the  best-fit  line  through  the  dry  data  is  almost 
horizontal.  The  data  at  60  percent  compaction  saturation  take 
the  form  described  above,  but  no  transition  point  can  be 
located  for  the  data  at  20  percent  compaction  saturation  (The 
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best-fit  line  for  the  data  at  2C  percent  compaction  saturation 
does  not  include  the  point  at  1  G,  the  lowest  gravity  scaled 
yield) . 

C.  CHATER  DIMENSIONS  VERSUS  DEPTH  OF  BURIAL 

The  graphs  of  apparent  crater  diameter  versus  depth  of 
burial  (DOB) ,  and  apparent  crater  depth  versus  DOB  are 
presented  in  Figures  8.9  and  8.10.  The  data  have  been 
normalized  to  remove  the  influence  of  charge  mass.  Data  from 
prototype  tests  at  Colorado  State  University  (CSU)  and  the 
full  scale  Sky  X  test  area^  at  Tyndall  AFB  are  included.  The 
tests  in  dry  sand  plot  below  the  tests  in  the  partially 
saturated  sand.  lie  trend  .issociated  with  the  degree  of 
saturation  during  compaction  is  apparent. 

The  best-fit  curve  for  sand  at  0  <  X  <  2 . 0  generated 
from  field  tests  at  WES  (Strange  et  al .  ,  1961)  is  also 
included  in  Figures  8.8  and  8.10.  These  curves  show  an 
optimum  DOE  in  the  range  of  1.25  to  1.5.  The  centrifuge  tests 
of  the  current  research  in  dry  and  partially  saturated  sand 
appear  to  indicate  an  optimum  DOB  for  apparent  crater  diameter 
and  depth  near  1.75. 

Figure  8.11  is  a  reproduction  of  Figure  V-10  of  the  Air 
Force  Protective  Construction  Design  Manual.  The  figure  shows 
the  predicted  range  in  apparent  crater  dimensions  from  dry 

^Personal  communication,  Walter  Buchholtz,  Structural  Research 
Engineer,  HQ  AFCESA/RACS,  Tyndall  AFB,  FL  32403. 
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Apparent  Crater  Diameter  (ft)  /  (Charge  Weight,  kg  TNT)l/^ 


Depth  of  Burial  (m)  /  (Charge  Mass,  kg  TNT)^^^ 


0.0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6 


Depth  of  Burial  (ft)  /  (Charge  Weight,  lb  TNT)^/'"^ 


Figure  8.9  Normalized  Apparent  Inside  Crater  Diameter 

Compared  to  the  U.  S.  Army  Engineer  'Waterways 
Experiment  Station  Prediction  Curve  (Strange, 
1961).  (Conversion:  0 . 3967f t/lb^/^  =  1.0 


Apparent  Croter  Diameter  (m)  /  (Charge  Moss,  kg  TNT)^/^ 


Apparent  Crater  Depth  (ft)  /  (Charge  Weight,  lb  TNT)l/3 


Depth  of  Burial  (m)  /  (Charge  Mass,  kg  TNT)!/^ 


0.0  0.2  0.4-  0.6  0.8  1.0  1.2  1.4  1.6 


Figure  8.10  Normalized  Apparent  Crater  Depth  Compared  to 
the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  Prediction  Curve  (Strange,  1961). 
(Conversion:  0.3967ft/lb^/^  =  l.Om/kgf/^) 


Apparent  Crater  Diameter  (m)  /  (Charge  Mass,  kg  TNT)l/^ 


Depth  of  Burial  (m)  /  (Charge  Moss,  kg  TNT)^/^ 
0.0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6 


Figure  8.11  Normalized  Apparent  Crater  Depth  and 

Inside  Diameter  Compared  to  the  .Air  Force 
Protective  Construction  Design  Manual 
Prediction  Curves. 

(Conversion:  0 . 39  67  f  t/lb^/ ^  =  1.0,m/kg^/^) 


Apparent  Crater  Diameter  (m)  /  (Charge  Mass,  kg  TNT)^/^ 


(S  ■-  0%)  to  wet  (S  >  0%)  sand  for  different  normalized  depths. 
The  Air  Force  design  curves  were  generated  using  the 
conventional  weapons  computer  program  (CONWEP,  17  February 
1988  version)  which  references  Figure  V-10. 

Two  observations  concerning  apparent  crater  diameter  are 
evident  from  Figure  8.11.  First,  the  majority  of  the  field 
and  centrifuge  data  lie  above  the  range  predicted  by  the 
CONWEP  program.  Second,  the  WES  (Strange  et  al.,  1961)  best- 
fit  curve  also  lies  above  the  range  predicted  by  the  CONWEP 
program.  The  WES  curve  correlates  better  with  centrifuge  and 
full  scale  data  than  the  curve  generated  by  the  CONWEP 
computer  program. 

In  contrast,  the  range  predicted  by  the  CONWEP  computer 
program  for  apparent  crater  depth  appears  to  fit  the  test  data 
well.  The  range  contains  the  central  data  from  the  Sky  X  and 
centrifug:  tests,  and  encompasses  the  WES  curve. 

Figure  8.12  shows  graphical  comparisons  of  the  current 
research  data  with  the  apparent  and  true  crater  dimensions  in 
"dry"  and  "wet"  sand  predicted  by  the  CONWEP  computer  program. 
As  in  Figure  8.11,  the  data  are  normalized  to  remove  the 
influence  of  charge  mass.  These  figures  illustrate  how  the 
apparent  crater  dimensions  of  the  centiifuge  tests  fit  the 
prediction  curves  of  the  Air  Force  Protective  Construction 
Design  Manual . 

The  scaled  apparent  crater  depth  data  trum  the  dry  tests 
of  the  current  research  follow  the  appaient.  crater  depth  curve 
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Crater  Dimensions  (m)/{Charge  Moss,  kg  TNT)''^'’  Crater  Dimensions  (m)/(Charge 


Charge  Boriol  (m)  /  (Charge  Moss,  kg  TNT; 

(b) 

Figure  8.12  Normalized  Apparent  Crater  Depcu  and  Inside 
Diameter  Compared  to  the  Apparent  and  True 
Crater  Dimensions  Predicted  by  the  Air  Force 
Conventional  Weapons  Computer  Program  (CONWEP) 
a)  Dry  Sand;  b)  Partially  Saturated  Sand. 
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predicted  for  "dry"  sand,  whereas  the  scaled  apparent  crater 
depth  data  from  the  partially  saturated  tests  vary  from  values 
less  than  the  apparent  depth  curve  to  values  greater  than  that 
true  depth  curve  predicted  for  "wet"  sand.  Due  to  the  unknown 
moisture  content  of  the  sand  data  used  by  the  CONWEP  computer 
program,  the  accuracy  of  the  partially  saturated  apparent 
crater  depth  data  can  not  be  assessed. 

The  majority  of  the  scaled  apparent  crater  diameter  data 
measured  in  the  current  research  are  larger  that  the  predicted 
true  crater  diameter  in  both  the  "dry"  and  "wet"  sand.  Thus, 
the  shape  of  the  craters  produced  on  the  centrifuge  is 
different  than  the  shape  of  the  craters  predicted  by  the 
CONWEP  coniputer  program  (Figure  8.13). 

A  best-fit  curve  from  the  DOB  analysis  by  Sclunidt  (1979) 
is  superimposed  on  the  data  from  the  current  research  (Figure 
8.14).  The  data  match  the  curve  at  the  lower  scaled  crater 
efficiencies  defined  as  (see  section  III)  . 
Unfortunately,  these  particular  data  points  have  the  largest 
confidence  intervals  (see  section  VII)  .  The  data  points 
having  the  smallest  confidence  intervals  are  the  points 
furthest  from  the  curve.  However,  data  collected  at  similar 
G-levels  (25  G  -  65  G)  by  Viktorov  and  Stepenov  (Schmidt, 
1979)  are  also  not  well  predicted  at  the  higher-scaled  crater 
efficiencies.  Thus,  the  data  from  the  current  research  are  in 
reasonable  agreement  with  the  best-fit  curve  derived  by 
Schmidt  (1979)  . 
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Figure  8.13  Crater  Profile  Comparison  Between  the 

Centrifuge  Apparent  Crater  and  the  Crater 
Generated  by  the  CONWEP  Computer  Program, 


D.  DESATURA.TION  ANALYSIS 

The  moisture  content  profile  in  the  sand  during  the  test 
is  analyzed  in  two  ways.  First,  the  desaturation  curve  as 
shown  in  Appendix  E  is  scaled  to  model  dimensions  (divided  by 
19  G  and  67  G) .  Second,  the  bucket  elevations  are  scaled  to 
field  elevations  (multiplied  by  19  G  and  67  G)  allowing 
comparison  to  desaturation  occurring  in  the  field. 

Using  the  first  approach,  the  moisture  content  analysis 
shows  that  desaturaticn  is  the  greatest  in  the  top  layer  of 
sand,  and  decreases  with  increasing  depth.  The  higher  the  G- 
level  attained  during  the  test,  the  greater  the  desaturation. 
At  20  percent  compaction  saturation  the  top  10  cm  of  sand 
desaturated  to  approximately  the  residual  water  content  of 
approximately  10  percent  at  all  G-levels  (Figure  8.15). 
Similar  trends  are  shown  for  30.5,  60,  and  70  percent 
compaction  saturation. 

The  second  approach  for  analyzing  the  moisture  content 
profile  allows  the  direct  comparison  of  moisture  content  data 
with  the  water  retention  curve  for  Tyndall  Beach  Sand.  In 
Figures  8.16  through  8.19,  which  show  compaction  saturation 
plotted  against  capillary  pressure,  the  centrifuge  data  have 
been  separated  by  G-level,  The  formation  of  a  water  table 
within  the  specimen  bucket  due  to  desaturation  during 
centrifugation  shifts  the  water  retention  curve  upward. 
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Saturation  (*) 


Capillary  Pressure  Versus  Saturation, 

19  G.  (Converion;  10.200171^20  “  l.OkPa) 
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Capillary  Pressure  (KPa) 


Figure  8.17  Capillary  Pressure  Versus  Saturation, 

26  G.  (Conversion:  10.20cmH2O  "  l.OkPa) 
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Copillary  Pressure  (KPo) 


Figure  8.18  Capillary  Pressure  Versus  Saturation, 

49  G.  (Conversion:  lO.ZOcniH^O  ”  l.OkPa) 
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E .  CRATER  DIMENSION  VERSUS  SATURATION  AT  COMPACTION 

In  Figure  8.20  the  scaled  apparent  crater  dimensions  are 
plotted  by  G-level  against  compaction  saturation.  As  stated 
previously  (section  V)  ,  the  tests  at  19  G  and  26  G  were 
designed  to  use  the  modeling  of  models  technique.  The  tests 
use  different  charge  sizes  to  model  craters  produced  by  7  kg 
of  TNT.  Thus,  the  apparent  crater  dimensions  from  these  tests 
should  be  identical.  If  the  apparent  crater  dimensions  are 
significantly  different,  then  a  charge  size  effect  exists. 

Figure  8 . 20  shows  the  apparent  crater  dimensions  for  the 
tests  at  19  G  and  26  G  are  the  same  within  the  bounds  of 
experimental  error.  As  Goodings  et  al.  (1938)  and  Schmidt  and 
Holsapple  (1978)  show,  no  charge  size  effect  exists  and  the 
modeling  of  models  is  valid. 

The  slight  rise  of  the  26  G  test  crater  radius  and  depth 
at  the  compaction  saturation  of  35  percent  (test  number 
120.63-40-1)  may  be  explained  by  experimental  procedure. 
Because  the  payload  platform  had  not  rotated  the  full  90 
degrees  at  26  G,  the  centrifuge  speed  was  increased  to  43  G  at 
which  point  full  rotation  occurred.  With  the  payload  platform 
held  in  the  vertical  position,  at  a  90  degree  rotation,  the 
centrifuge  was  slowed  to  26  G. 

Figure  8. 15. a  shows  that  for  tests  at  40  percent 
compaction  saturation,  the  saturation  at  the  detonator 
elevation  at  the  conclusion  of  the  test  was  significantly 
lower  at  26  G  (S  =  9  percent,  long  dash)  than  at  19  G  and  49 
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Scoled  Apporenl  Croler  Depth  (m)  Scoled  Appdreot  Croter  f-^odius  (m)  Scaled  Apporent  Croter  Volume  (m^) 


(a) 


(b) 


Figure  8.20  Ar^arent  Crater  Dimensions  Versus  Ccmpaction 
Siiuration:  a)  Scaled  Dimensions; 
o’  Normalized  Dimensions. 
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G  (S  =  12  percent,  solid  line  and  medium  dash  respectively). 
In  moist  sand,  a  lower  pore-water  pressure  induced  by  the 
lower  saturation  increases  the  apparent  cohesion  component  of 
shear  strength.  The  equation  given  by  Fredlund  (1986)  for  the 
shear  strength  of  unsaturated  soil  is 

s  =  -  u^)  tan4>*  +  (o„  -  u^)  tari(j>' 

=  Apparent  Cohesion 

^  Frictional  Shear  Resistance  (8.1) 

where  u,  is  the  pore-air  pressure,  u,.  is  the  pore-water 
pressure,  (J)*^  is  the  angle  of  shear  strength  increase  with  an 
increase  in  (u,  -  u^)  ,  is  the  total  effective  normal  stress, 
and  <{>'  is  the  effective  angle  of  internal  friction. 

For  the  centrifuge  tests,  u,  was  at  atmospheric  pressure 
(u,  =  0)  .  In  unsaturated  soils,  u^  is  always  less  than  the  air 
pressure,  u,,  thus  u,  must  be  negative  (Fredlund,  1986).  For 
dry  soil,  u,  is  zero,  becoming  negative  as  moisture  is  added 
to  the  soil  (Fredlund,  1986) .  However,  as  the  soil  moisture 
increases,  u^  becomes  less  negative.  Therefore, 
theoretically,  the  apparent  crater  dimen.sions  of  the  26  G 
specimen  spun  to  43  G  should  be  smaller  than  the  apparent 
crater  dimensions  of  the  19  G  specimen  because  the  sand  had  a 
lower  saturation  and  therefore  a  higher  shear  strength.  An 
equivalent  explosive  mass  should  have  ejected  more  sand  from 
the  weaker  19  G  specimen  than  the  stronger  26  G  specimen. 
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The  test  at  a  compaction  saturation  of  20  percent  and  67 
G  (test  number  191.89-20-40)  was  also  spun  to  a  higher  G-level 
(9.1  G)  to  correct  the  orientation  of  the  payload  platform. 
However,  in  this  case,  the  apparent  crater  dimensions  at  67  G 
increase  and  decrease  at  the  same  compaction  saturation  levels 
as  the  apparent  crater  dimensions  at  49  G.  By  67  G,  the  sand 
at  the  detonator  had  already  been  desaturated  to  the  resiJual 
level  (9  to  30  percent). 

Figure  8.20  also  shows  that  scaled  apparent  crater 
volume,  radius,  and  depth  increase  with  G-level.  When  the  G- 
level  is  multiplied  by  a  factor  of  2.5  (from  26  G  to  67  G) 
which  corresponds  to  increasing  the  prototype  explosive  charge 
from  7  to  124  kg,  apparent  crater  volume  increases  by  an 
order  of  magnitude  from  7  m-  to  70  m--  and  doubles  apparent 
crater  radius  from  2,3  ro  to  4.6  m. 

The  scaled  aj)parent  crater  volume,  radius,  and  depth 
remain  constant  from  20  to  70  percent  compaction  saturation. 
A  noticeable  increase  in  scaled  apparent  crater  volume  and 
depth  is  observed  between  0  and  20  percent  compaction 
saturation.  An  increase  in  scaled  apparent  crater  dimensions 
is  expected  to  occur  a.*:  low  as  the  residual  level  of 
saturation  where  the  capillary  forces  are  the  largest. 

An  increase  in  stress  transmission  has  been  observed  to 
occur  in  Ottawa  Sand  between  0  and  the  residual  level  of 
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compaction  saturation.^  As  shown  by  Fredlund  (1986),  the 
apparent  cohesion  component  of  soil  shear  strength  increases 
as  the  pore-water  pressure  decreases,  becomes  more  negative. 
Thus,  in  that  cratering  is  controlled  by  soil  strength,  the 
effects  of  moisture  should  increase  as  the  degree  of 
saturation  approaches  the  residual  level.  Conversely,  at 
saturations  greater  than  the  residual  level,  the  less  effect 
capillarity  should  have  on  crater  formation. 

F.  ASPECT  RATIO  VERSUS  SATURATION  AT  COMPACTION 

To  analyze  crater  shape,  the  aspect  ratio,  as  defined  as 
apparent  crater  radius  divided  by  apparent  crater  depth  (r /d) , 
was  plotted  against  compaction  saturation  (Figure  8.21) .  The 
data  are  separated  by  G-level.  The  observations  from  this 
graph  are;  1)  craters  in  dry  sand  have  a  larger  aspect  ratio 
at  all  G-levels;  2)  charge  size  does  not  alter  aspect  ratio; 
and  3)  the  addition  of  moisture  significantly  reduces  the 
aspect  ratio.  Also,  based  on  a  few  tests  at  1602  )cg/m^,  there 
is  strong  indication  that  under  dry  conditions,  an  increase  in 
density  will  lower  the  aspect  ratio. 

Figure  8.21  reinforces  the  influence  of  water  in  the 
specimen.  The  average  aspect  ratio  of  the  dry  specimens  is 
3.17  with  a  range  between  2.82  and  3.40,  whereas  the  aspect 
ratios  of  the  partially  saturated  specimens  range  from  1.3  to 

^Personal  Communication,  C.  Allen  Ross,  Professor,  University 
of  Florida,  Gainsville,  Florida;  Consultant,  HQ  AFCESA/RAC.S, 
Tyndall  AFB,  Florida. 
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Aspect  Ratio,  r/d 


Figure  8.21  Aspect  Ratio  Versus  Compaction  Saturation. 


2.0.  Lines  of  average  aspect  ratios  are  drawn  through  the 
partially  saturated  data.  Table  8.5  lists  the  average  aspect 
ratio  and  the  apparent  crater  wall  angle  for  each  G-level . 
The  crater  wall  angle  (cx)  is  calculated  from  the  equation  tan 
a  =  d/r  (Figure  8.22) . 


TABLE  8.5.  AVERAGE  ASPECT  RATIO  AND  APPARENT  CRATER  WALL 
ANGLE  AS  TAKEN  FROM  FIGURE  8.21. 


S. 

C 

9 

PBX 

(gm) 

Scaled 
TNT  EQ 
(KQ) 

.  DOB 
(cm) 

r/d 

a 

(degl 

17 

to 

70 

18. 

86 

1 

.043 

7. 

6 

7.6 

1.3 

37.6 

17 

to 

70 

26. 

34 

0 

.383 

7. 

6 

5.44 

1.4 

35.5 

20 

to 

60 

49. 

.11 

0 

.383 

49. 

4 

5.44 

1 . 9 

27.8 

20 

to 

60 

66. 

,65 

0 

.383 

123. 

6 

5.44 

2.0 

26.6 

0 

18. 

,86 

1 

.043 

7  . 

6 

7.6 

3.2 

17.5 

26, 

,34 

0 

.383 

7 . 

6 

5.44 

Ave . 

Ave . 

49. 

,11 

0 

.383 

49. 

4 

5 . 44 

66. 

.65 

0 

.383 

123. 

6 

5.44 

One  possible  explanation  for  the  high  aspect  ratio  for 
the  craters  in  dry  sand  is  slope  failure.  Slumping  of  the 
side  wall  would  decrease  the  crater  depth,  increase  the  crater 
radius,  and  raise  the  aspect  ratio.  The  measured  apparent 
crater  depth  is  shallower  in  dry  craters  (Figure  8.23),  and 
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Figure  8.23  Measured  Apparent  Crater  Depth  Versus 
Saturation.  All  Data  are  from  Charge 
an  Equivalent  of  383  mg  of  PBX . 


the  crater  wall  angle  of  17.5®  (Table  8.5)  is  much  lower 
thanthe  angle  of  repose  for  beach  sand  (approximately  30°-35°)  . 

With  moisture,  sand  develops  capillary  forces  which 
create  an  apparent  cohesion.  Evidence  in  support  of  the 
capillary  forces  and  an  apparent  cohesion  is  provided  by  the 
high  crater  wall  angles  of  37.6®  and  35.5®  calculated  for  the 
specimens  at  19  G  and  26  G  (Table  8.5).  Higher  wall  angles 
permit  deeper  craters  which  result  in  lower  aspect  ratios. 
Thus,  slope  failure  did  not  occur  in  the  craters  formed  at  19 
G  and  26  G. 

Figure  8.20  shows  apparent  crater  depth  increases  more 
between  degrees  of  compaction  saturation  of  O^and  20  percent 
than  in  apparent  crater  radius.  As  previously  discussed,  20 
percent  compaction  saturation  is  not  expected  to  be  the 
limiting  degree  of  compaction  saturation.  Due  to  the  rise  in 
capillarity  with  decreasing  moisture,  the  residual  level  of 
saturation  is  instead  expected  to  be  the  limit  to  increasing 
apparent  crater  dimensions. 

As  the  G-level  increases  to  49  G  and  67  G,  the  apparent 
cohesion  component  of  sand  shear  strength  becomes  a  smaller 
percentage  of  the  total  strength  (Figure  8.24),  and  slope 
failure  reoccurs.  The  measured  apparent  crater  depth  for 
tests  at  49  G  and  66  G  is  less  than  the  DOB  of  5.44  cm  (Figure 
8.23),  and  the  crater  wall  angles  of  27.8°  and  26.6°  (Table 
8.5)  are  significantly  lower  than  the  crater  wall  angles  for 
tests  at  19  G  and  26  G.  For  degrees  of  compaction  saturation 
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Cohesion  =  1  0G/{ 


Figure  8.24  Decrease  of  the  Cohesion  Component  of  Soil 
Strength  Relative  to  the  Total  Shear 
Strength  with  an  Increase  in  G-Level 
(Stress ) . 
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from  20  to  60  percent,  the  apparent  crater  depth  increases 
less  than  the  apparent  crater  radius  (Figure  8.20)  which  leads 
to  the  greater  aspect  ratios  shown  in  Figure  8.21  for  the 
tests  at  49  G  and  67  G.  Figure  8.25  depicts  these  differences 
in  aspect  ratio  for  the  tests  at  26  G,  49  G,  and  67  G  through 
the  crater  profile. 

Figures  8.26  and  8.27  show  how  subjecting  specimens  to 
elevated  G' s  prior  to  a  test  at  1  G  raises  the  aspect  latio. 
The  apparent  craters  become  proportionally  shallower,  and  both 
the  apparent  crater  depth  and  apparent  crater  diameter 
decrease  (Table  8.6).  Moisture  is  drained  from  the  detonator 
elevation  which  increases  the  capillary  forces  and  in  turn 
increases  the  soil  strength. 
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Aspect  Rotio,  r/d 


2.5 


2.0 


Compoction  Saturation  (??) 


Figure  8.26  Differences  in  Aspect  Ratio  due  to 

Centrifugation  Prior  to  Testing  at  1  G  . 
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TABLE  8.6.  THE  DECREASE  IN  ACTUAL  APPARENT  CRATER 

DIMENSIONS  DUE  TO  CENTRIFUGATION  AT  57  G 
(179.5  RPM)  PRIOR  TO  TESTING  AT  1  G. 


Inside 


Test  Number 
(rpm-lG (sat) 

-#) 

Diameter 

(cm) 

Depth 

(cm) 

Volume 

(cm^) 

0.0-lG(20) 

-1 

24 . 13 

7.78 

1368 

179.5-1G(20) 

-3 

23.02 

6.48 

1038 

0.0-30(60) 

-2 

25 . 40 

9.17 

1869 

179. 5-lG(60) 

-A 

24 . 45 

7 . 94 

1427 

G.  CORIOLIS  EFFECT  ON  CRATER  SHAPE 

The  asymmetry  of  the  crater  molds  pictured  in  Appendix  A 
is  evidence  that  the  use  of  centripetal  acceleration  distorts 
the  shape  of  the  apparent  crater  in  sand.  In  every  instance, 
the  distortion  occurs  towards  the  We.st.  Figure  8.28  shows  the 
zone  in  which  sand  grains  have  been  broken  into  fine  particles 
by  the  blast  is  centered  about  the  bottom  of  the  apparent 
crater.  Thus,  there  is  an  indication  the  true  crater  as  well 
as  the  apparent  crater  is  shifted  towards  the  west 

In  field  tests,  crater  skewness  is  usually  attributed  to 
a  tilted  charge;  a  charge  whose  longitudinal  axis  is  not 
vertical  at  the  time  of  detonation.  However,  the  bottom  of 
all  the  craters  formed  on  the  centrifuge  have  been  visibly 
skewed  in  the  West  direction,  opposite  the  rotation  of  the 
centrifuge.  Figure  8.29  shows  the  profile  of  test  number 
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Figure  8.28  Typical  Location  of  Crushed  Sand  Grains 
in  the  East -West  Direction: 
a)  Photograph;  b)  Schematic. 


Figure  8.29  Crater  Mold  Prefixes  in  the  Uast-West  and 
North-South  Directions. 
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191.89-60-3  taken  in  the  East-West  and  North-South  directions. 
The  procedure  used  to  generate  the  profile  is  described  in 
Appendix  G.  In  the  North-South  direction,  the  apparent  crater 
is  symmetric  about  the  axis  of  the  detonator.  If  the  Coriolis 
acceleration  were  not  the  cause  of  the  distortion  of  the 
apparent  crater  shape,  the  asymmetry  would  not  be  confined  to 
the  East— West  direction.  Thus,  the  symmetry  of  the  apparent 
crater  in  the  North-South  direction  is  further  verification 
that  the  Coriolis  acceleration  is  altering  the  shape  of  the 
apparent  crater.  The  angle  between  the  lines  drawn  from  the 
intersection  of  the  plane  containing  the  detonator  and  the 
plane  of  the  ground  surface  to  the  lowest  point  of  each  crater 
profile  is  15.5°. 

The  skewness  of  several  other  craters  was  visually 
approximated.  These  approximations  are  listed  in  Table  8.7. 
The  data  indicate  that  greater  G-levels  produce  greater 
apparent  crater  asymmetry  (Appendix  I) ,  however,  additional 
research  is  required. 

Schmidt  et  al .  (1986)  and  Steedman  (1990)  attribute 
apparent  crater  asymmetry  to  the  uneven  distribution  of 
ejecta,  but  the  Coriolis  acceleration  should  effect  more  than 
just  the  ejecta.  Because  of  high  initial  velocities  produced 
by  the  blast,  soil  particles  moving  in  directions  other  than 
the  North-South  direction  should  experience  a  Coriolis  force. 
Thus,  a  vectorial  analysis  was  performed  to  show  how  the 
Coriolis  acceleration  (a^)  effects  the  rupture  zone  as  well  as 
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TABLE  8.7.  APPARENT  CRATER  ASYMMETRY  DUE  TO  THE  CORIOLIS 
ACCELERATION. 


Test 

(rom-sat-#) 

G-level 

(G) 

Angle  of  Rotation 
{±2  dec) 

1C2.0C  70-3 

18. 8C 

11.0 

120.63-70-4 

26.34 

13.0 

120.63-60-1 

26.34 

12.8 

164.71-20-5 

49.11 

15.7 

164.71-40-4 

49.11 

19.8 

164.71-60-4 

49.11 

20.0 

191.89-20-4 

66.65 

18.6 

191.89-40-3 

66.65 

16.0 

191.89-60-3 

66.65 

15.5 

the  ejecta. 

The  blast  is  assumed  to  be  a  point  source  and  the 
velocity  vectors  of  the  particles  being  ejected  from  the 
centrifuge  bucket  4re  separated  into  radial  and  tangential 
components.  The  Coriolis  acceleration  vectors  are  the  vector 
cross  products  between  the  centrifuge  rotation  vector  and  the 
relative  particle  velocity  vectors.  The  Coriolis  force 
vector,  however,  will  act  in  the  direction  opposite  that  of 
the  Coriolis  acceleration.  The  inertia  force  opposing  the 
Coriolis  acceleration  will  govern  the  movement  of  the  soil 
particles . 

Figure  8.30  shows  the  directions  of  the  Coriolis  force 
components  drawn  in  the  plane  of  the  original  ground  surface 
about  the  edge  of  a  spherical  crater.  The  Coriolis  force 
vectors  derived  from  the  tangential  velocity  of  the  particles 
point  out  of  the  bucket  on  the  West  side  of  the  crater  and 
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into  the  bucket  on  the  East  side  of  the  crater.  All  of  the 
Coriolis  force  vectors  derived  from  the  radial  velocity  of  the 
particles  point  to  the  East.  The  Coriolis  force  acting  on  the 
soil  below  the  detonator  will  be  pointed  to  the  West,  assiiming 
the  Coriolis  effect  is  reversed  for  particles  accelerating 
down  into  the  bucket.  Therefore,  during  the  explosion,  the 
Coriolis  forces  tend  to  rotate  the  soil  mass  about  the  North- 
South  axis  at  the  detonator.  The  rotation  of  the  soil  mass 
surrounding  the  detonator  by  the  Coriolis  acceleration  will 
shift  the  bottom  of  the  crater  to  the  West  against  the 
centrifuge  rotation.  Thus,  the  analytical  result  supports  the 
skewed  crater  shape  produced  on  the  centrifuge. 


127 


SECTION  IX 


CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

1 .  SCALING  LAWS  FOR  DRY  AND  MOIST  SAND 

Craters  induced  by  explosive  detonations  in  dry  and 
partially  saturated  Tyndall  Beach  Sand  were  modeled  in  a 
centrifuge.  The  apparent  dimensions  of  the  explosion-induced 
craters  are  analyzed  as  dimensional  terms  and  dimensionless  pi 
terms . 

In  the  analysis  of  dimensional  terms,  the  scaled  apparent 
crater  volume,  radius,  and  depth  are  plotted  against  the  TNT 
equivalent  of  the  scaled  explosive  mass  (Figures  8.1  to  8.6). 
The  analysis  shows  through  a  comparison  of  yield  exponents 
that  the  scaling  laws  for  apparent  crater  dimensions  are  the 
same  for  dry  and  partially  saturated  sand  (Table  9.1) .  The 
yield  exponents  of  the  current  research  are  also  similai  to 
the  yield  exponents  of  crater  dimensions  produced  on 
centrifuges  by  Goodings  et  al .  (1988)  and  Schmidt  and 
Holsapple  (1978,  1980),  and  produced  in  the  field  by  the  U.S. 
Army  Engineer  Waterways  Experiment  Station  (WES)  (Strange  et 
al.,  1961;  Rooke  et  al . ,  1974). 
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Radius  -  Moist  0.308  0-316 _ 0.306 _ 0-305 _ 

Depth  -  Dry  0.279  0.28 _ 0.279 _ 0.2'10 

Depth  -  Moist  0.297  0-335 _ 0.238 _ 

-  Yield  exponents  for  the  dry  test  data  do  not  include  1  G  test  results. 

Yield  exponents  for  the  moist  test  data  include  1  G  test  results. 


Conversely,  in  the  analysis  of  the  pi  terms  derived  by 
Schmidt  and  Holsapple  (1977)  for  buried  charges,  the  apparent 
crater  volume  (K,  =  vp/W)  and  depth  (n<j  =  d(p/W)^'^)  data  at 
compaction  saturations  of  20  and  60  percent  do  not  parallel 
the  data  for  dry  beach  sand  when  plotted  in  Figure  8 . 7  against 
the  gravity  scaled  yield  (rtj  =  (G/Q)  (W/5)  .  The  apparent 
crater  volume  and  depth  data  support  the  material  strength 
model  developed  by  Schmidt  and  Holsapple  (1979).  At  low 
gravity  scaled  yields  the  volume  and  depth  curves  are 
horizontal;  the  apparent  cohesion  of  the  particlly  saturated 
sand  dominates  the  cratering  process.  The  lower  the 
compaction  saturation  and  the  larger  the  apparent  cohesion, 
the  higher  the  transition  point  becomes  between  cohesion 
dominated  cratering  and  frictional  shear  resistance  dominated 
cratering. 

Thus,  when  dimensional  terms  are  analyzed,  the  scaled 
apparent  crater  dimensions  of  volume,  radius,  and  depth  follow 
the  same  scaling  relations  in  partially  saturated  sand  as  in 
dry  sand.  However,  when  the  dimensionless  pi  terms  are 
analyzed,  the  apparent  crater  volume  and  depth  do  not  follow 
the  same  scaling  relations  in  partially  saturated  sand  as  in 
dry  sand.  No  conclusions  can  be  drawn  from  the  pi  term  for 
apparent  crater  radius. 
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2 .  DEPTH  OF  BURIAL  EFFECT 

In  a  second  analysis  of  dimensional  terms,  the 
apparent  crater  dimensions  are  normalized  (divided  by  the  cube 
root  of  the  explosive  mass)  and  plotted  against  the  normalized 
depth  of  charge  burial  (DOB) .  The  normalized  apparent  crater 
diameter  and  depth  data  follow  the  best-fit  prediction  curves 
given  in  Figures  8.9  and  8.10  for  sand  derived  by  WES  (Strange 
et  al.,  1961) .  The  craters  dimensions  measured  in  dry  sand 
plot  below  the  crater  dimensions  measured  in  partially 
saturated  sand.  Amongst  the  partially  saturated  data,  no 
trends  due  to  the  differences  in  compaction  saturation  are 
observed . 

The  apparent  crater  diameter  data  do  not  fall  within  the 
range  from  dry  (field  dry)  to  wet  (saturated)  sand  predicted 
by  the  Air  Force  Protective  Construction  Design  Manual  and  the 
Army  design  manual,  Fundamentals  of  Protective  Design  for 
Conventional  Weapons  (Figure  8.11).  The  scaled  apparent 
crater  diameter  data  for  craters  formed  in  dry  and  partially 
saturated  Tyndall  Beach  Sand  are  approximately  50  percent 
greater  than  the  predicted  apparent  and  true  crater  dimensions 
in  "dry"  and  "wet"  sand  (Figure  8.12). 

The  scaled  apparent  crater  depth  data  for  craters  formed 
in  dry  sand  follow  the  curve  predicted  by  the  Air  Force 
Protective  Construction  Design  Manual  and  the  Army  design 
manual,  Fundamentals  of  Protective  Design  for  Conventional 
Weapons,  for  apparent  crater  depth  in  "dry"  sand  (Figure 
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8. 12. a).  In  "wet"  sand,  the  scaled  apparent  crater  depth  data 
at  a  normalized  DOB  of  0.73  m/kg^^^  bracket  the  predicted  range 
for  apparent  and  true  crater  depth  (Figure  8.12.b). 

3 .  DESATURATION  DURING  CENTRIFUGATION 
Centrifugation  is  shown  to  cause  desat ’"^ation  in 

sand  ranging  in  degrees  of  saturation  from  20  percent 
(Figures  8.18).  The  desaturation  is  such  thac  the  water 
retention  curve  is  reproduced  within  the  specimen  bucket 
(Appendix  E) .  The  greatest  desaturation  occurs  at  the  top  of 
the  sample,  and  the  higher  the  g-level,  the  greater  the 
desaturation.  The  results  indicate  that  capillary  rise  scales 
to  the  prototype  saturation  profile  as  1/N,  where  N  is  tie 
model  scale. 

4 .  SATURATION  VERSUS  CRATER  DIMENSION 

The  scaled  apparent  crater  dimensions  of  crater 
volume,  radius,  and  depth  increase  by  150,  20,  and  100  percent 
respectively  between  0  and  20  percent  compaction  saturation 
(Figures  8.20  and  8.23).  If  the  influence  of  moisture  on 
stress  transmission  ratios  can  be  applied  to  cratering,  the 
increase  in  scaled  apparent  crater  dimensions  will  occur  at 
the  residual  level  of  saturation.  In  Tyndall  Beach  Sand,  the 
residual  level  of  saturation  is  approximately  10  percent. 
Beyond  20  percent,  the  compaction  saturation  is  not  shown  to 
influence  apparent  crater  dimensions. 
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5 .  EFFECT  OF  MOISTURE  ON  APPARENT  CRATER  ASPECT  RATIO 
Craters  in  dry  Tyndall  Beach  Sand  have  aspect  ratios 
of  approximately  3.2.  An  increase  in  the  unit  weight  of  the 
material  or  the  addition  of  moisture  lowers  the  aspect  ratio. 
At  the  residual  level  of  saturation,  the  aspect  ratio  falls  to 
1.3  and  1.4  for  19  G  and  2  6  G,  respectively.  At  the  higher  G- 
levels  of  49  G  and  67  G,  the  aspect  ratio  increases  to  2.0 
(Figure  8.21)  . 

The  crater  wall  angles  of  the  craters  measured  in 
partially  saturated  sand  are  37.6°  and  35.5°  from  the 
horizontal  for  19  G  and  26  G,  respectively.  The  decrease  in 
apparent  crater  depth  for  tests  performed  at  49  G  and  67  G, 
and  the  relatively  low  crater  wall  angles  of  the  craters 
measured  in  dry  sand  (17.5°)  and  in  the  partially  saturated 
sand  at  49  G  and  67  G  (27.8°,  26.6°)  suggest  slope  failure  as 
the  mechanism  causing  large  aspect  ratios  for  craters  formed 
in  dry  sand  and  at  49  G  and  67  G.  Centrifugation  of  a 
specimen  prior  to  testing  at  1  G  is  also  shown  to  raise  the 
aspect  ratio. 

6.  CORIOLIS  EFFECT  ON  CRATER  SHAPE 

Craters  produced  on  a  centrifuge  by  a  buried  charge 
are  skewed  in  the  plane  of  rotation  against  the  direction  of 
rotation  (Figure  8.29).  The  skew  of  the  crater  is  attributed 
to  the  Coriolis  acceleration  on  particle  motion.  The  effect 
of  the  Coriolis  acceleration  on  the  entire  crater  is  to  rotate 
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the  soil  mass  of  the  crater  about  an  axis  centered  at  the 
charge  and  parallel  to  the  angular  rotation  vector  of  the 
centrifuge . 


B.  RECOMMENDATIONS 

Cratering  depends  on  soil  strength  and  in  partially 
saturated  sand,  the  apparent  cohesion  is  often  a  significant 
component  of  soil  strength.  Thus,  further  study  of  cratering 
in  partially  saturated  sand  should  begin  with  calculations  of 
the  apparent  cohesion.  An  apparent  cohesion/saturation  curve 
would  prove  useful  in  subscale  laboratory  tests  and  full  scale 
field  tests.  A  known  apparent  cohesion  would  permit  a  measure 
of  the  true  soil  strength  to  be  determined  and  the  graphical 
analysis  of  a  cohesion  pi  term. 

Secondly,  colored  sand  layers  should  be  placed  above  the 
detonator  to  map  the  distribution  of  ejecta.  The  relationship 
between  apparent  and  true  crater  dimensions  are  not  known  for 
centrifuge  crater  models.  The  relationship  may  or  may  not  be 
the  same  for  full  scale  field  tests. 

Lastly,  centrifuge  results  may  not  be  consistent 
throughout  all  G-levels.  The  data  in  Figures  8.9,  8.10  and 
8.11  indicate  that,  for  a  constant  DOB,  normalized  apparent 
crater  dimensions  decline  with  an  increase  in  G-level .  Data 
collected  by  Goodings  et  al .  (1988)  and  Schmidt  and  Holsapple 
(1980)  support  this  observation  (see  Appendix  L) .  Thus,  the 
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effects  of  centripetal  acceleration  for  explosion-induced 
crater  formation  requires  further  definition. 
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APPENDIX  A 

MEASURED,  SCALED,  AND  NORMALIZED  APPARENT  CRATER  DATA 


141 


CONSTANT  NORMALIZED  DOB  SERIES 


TEST  DATA  NOTES  FOR  TYNDALL  BEACH  SAND  COMPACTED  TO  A  DRY 
DENSITY  OF  95  Ib/ft^  (1521.7  kg/m^)  . 

Saturation  Note 

17%,  20%  Test  numbers  191.89-20-4  subjected  to  93-g 

prior  to  being  tested  at  67-g. 

40%  Test  numbers  102.08-40-1  and  102.08-40-2 

tested  at  a  saturation  of  35%,  not  a 
saturation  of  40%. 

Test  number  120.63-40-1  subjected  to  43-g 
prior  to  testing  at  26-g. 

70%  Test  102.08-70-1  had  1.0  in  of  water  in  the 

crater  immediately  after  it  was  formed. 
Subtracting  1.0  in  from  a  measured  crater 
depth  of  3.22  in  places  the  water  2.22  in 
(5.64  cm)  from  the  soil  surface.  Subtracting 
5.64  cm  from  the  total  soil  depth  of  28  cm 
left  a  elevation  of  22.4  cm.  Scaling  this 
elevation  a  prototype  dimension  requires 
multiplying  the  22.4  cm  by  (18.86)^  which  is 
421.73  cm. 

Test  120.63-70-2  had  0.55  in  of  water  in  the 
crater  immediately  after  it  was  formed. 
Subtracting  0.55  in  from  a  measured  crater 
depth  of  2.28  in  places  the  water  1.73  in 
(4.39  cm)  from  the  soil  surface.  Subtracting 
4.39  cm  from  the  total  soil  depth  of  26  cm 
left  a  elevation  of  21.6  cm.  Scaling  this 
elevation  a  prototype  dimension  requires 
multiplying  the  21.6  cm  by  (26.34)^  which  is 
569.1  cm . 

Test  102.08-70-3  and  120.83-70-4  prepared  with 
one  layer,  the  bottom  layer,  of  dry  sand  to 
prevent  a  rise  in  the  water  table. 


A  separate  test  series  in  Ottawa  sand  is  also  included.  The 
OT  in  the  test  number  indicates  Ottawa  sand.  In  this  series, 
the  dry  density  of  specimens  were  100.67  Ib/ft'  in  order  to 
keep  the  relative  density  equal  to  that  of  the  Tyndall  sand 
specimens  (39% )  . 
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MEASURED  APPARENT  CRATER  DIMENSIONS 
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Note:  Test  number  102.08-0-7  indicates  explosive  detonation  at  102.08  rpm, 

zeio  percent  saturation  series,  test  seven.  A  separate  test  series  in 
Ottawa  sand  is  also  listed.  The  OT  in  the  test  number  indicates  Ottawa 
sand.  *The  dry  density  of  the  Ottawa  sand  specimen  was  1601.8  kg/m^ 


SCALED  APPARENT  CRATER  DIMENSIONS 
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Dry  density  of  the  Ottawa  sand  specimen  was  1601.8  kg/m"^  (100  Ib/ft 


NORMALIZED  APPARENT  CRATER  DIMENSIONS 
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Dry  density  of  the  Ottawa  specimen  was  1601.8  kg/m  ( 100  Ib/ft 


4 


Specimen  spun  to  93  G  prior  to  test  at  66.65 
Dry  density  of  the  Ottawa  sand  specimen  was 
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Dry  density  of  the  Ottawa  sand  specimens  were  1601.8  kg/m  (100  Ib/ft 


Cfi 

U  Q 

H  2 


^  u 


cl  ^  i  ?  a:  s  ? 

^11  m  «  g  ^ 


t  -' 


go  ■©  ^ 

d  o  o  o  c> 


3^0  ^  m  ^ 

^  w  O  M 


^  I  -P 


o  Fa  N  o 
r4  r4  r4  f4 


ri  ♦  tt  ^  rn  ^ 


a  s  s  s  s  s 

K  r»  r-  r^  r* 

d  d  d  d  d  d 


^  ^  >0  >0 
5  a  s  a 


r*.  <N  'n  ^  ^  ^ 


r<Q^  —  r4  — 
O  ^  >6  ^  'O  ^ 


3  ^  S  ^ 

I  L  < 

^  w  ft 


s  ^  s  5?  CO 

«  I  I  5 


- 1  88«5Sg 

I  (Sl  ri  ^  ^  ^  ^ 


O'  q 
r4  r4  fn 


Koo  »n  r^  op  o 

r^  r^  ^  ^  M 

wS  ^r»  vS  wS  wS  d 


p  S!  a  s  a  2 

d  d  d  d  d  •>« 


^  sr>  vn  ^  »n 

m  m  f^.  m 

op  00  oo  op  op  00 


—  Ci  — •  o  :fi  o 

00  oo  o  W-)  do 

d  d  d  pi  ^  pi 

^  ra  ra  pi 


eg  m  ^  f5  r5 

a  aa  a  cc 


cr  ^  O'  oo  S' 
■o  oc  ^  3c 

^  ^ 


148 


Dry  densiti'  of  the  Ottawa  sand  specimens  were  1601.8  kg/m 

( 100  Ib/ft^)  , 
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Specimen  spun  to  43  G  prior  to  test  at  26.34 
Dry  density  of  the  Ottawa  sand  specimen  was 


SCALED  APPARENT  CRATER  DIMENSIONS 
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Dry  density  of  the  Ottawa  sand  specimen  was  1601.8  kg/m  (100  Ib/f 


NORMALIZED  APPARENT  CRATER  DIMENSIONS 
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MEASURED  APPARENT  CRATER  DIMENSIONS 
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SCALED  APPARENT  CRATER  DIMENSIONS 
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Dry  density  of  the  Ottawa  sand  specimen  was  1601,8  kg/m  (100  Ib/ft 


NORMALIZED  APPARENT  CRAfER  DIMENSIONS 
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Dry  density  of  the  Ottawa  sand  specimen  was  1601.8  kg/m  (100  Ib/ft 
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Test  numbers  3  &  4  were  prepared  using  2  inches  of  dry  sand  for  the  first 
bottom  -  layer. 


NORMALIZED  APPARENT  CRATER  DIMENSIONS 
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120.63  -  70  -  2  16.82  1.835  12.74  5.61  2.80  1.95  134.55 


VARIATION  IN  DOB  SERIES 
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201,76 -VO -5  0.00133  72.18  2.70  9  66  5.75  0.86  8.0  26.72 

139.30 -VO -6  0.00230  35.17  4.08  1169  7  88  0.98  16.8  9.82 

173.75  -VO -7  0.00133  57.29  3.40  9.41  5.50  0.61  13  34.91 

189.94  -VO -8  0.00181  65.07  3.77  10.69  6.63  0.67  7.9  73.13 
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201.76  -VO -5  500.00  545.00  162  4  58.08  34.59  17.29  5.17  1739.71 

1.39.30  - VO -6  100.00  109.00  11.96  34.26  23.08  11.54  288  42233 

178  75 -VO -7  250  00  27250  16.23  44  91  26.26  13  13  291  798  84 

189.94 -VO -8  500.00  545.00  20.44  57.%  35.92  17.%  3.64  1263.88 


NORMALIZED  APPARENT  CRATER  DIMENSIONS 

TYNDALL  BEACH  SAND 
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13930 -VO -6  109.00  Z503  7.17  4  83  2  42  0.60  88.41 

178.75  -VO -7  272.50  Z504  6.93  4.05  Z03  0.45  123.22 

18994  -VO -8  545  00  2  503  7.10  4  40  220  0.45  154.73 


1  G  SERIES 
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vibrated  - 1;  0  -  5  0.00133  1.00  0  2.49  10.50  625  0.8S  202  3329 
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vibrated  -  IjO- 5  0.0013  0.0014  0.2073  0.8750  03208  02604  0.0729  0.01167 

pluviated  -  lgO  -6  0.0013  0  0014  02073  08542  03000  02500  0.0573  0.00468 


NOR\IALIZED  APPARENT  CR.4TER  DIMENSIONS 


vibrated  IgC  - 5  0.30M  1A32  7.732  4.603  2  301  0  644  0.10311 

(>lu%isled-l{0  -  6  0.0014  l.a32  7.346  4.419  2209  0306  0.04132 
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MOISTURE  CONTENT-SATURATION  CONVERSION  CURVE 
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Saturation  (9S) 


Figure  C.l  Moisture  Content-Saturation  Conversion 
Curve  for  Tyndall  Beach  Sand. 
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MOISTURE  CONTENT  DATA 
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NOTES  FOR  MOISTURE  CONTENT  DATA. 


Test  No. 


Note 


120.63-40-1  Run  up  to  43-g  prior  to  testing  at  26-g 

(Balance  light  did  not  come  on) . 

One  side  of  the  bucket  consistently  wetter 
than  the  other,  12  o'clock  versus  3  o'clock. 
Moisture  contents  were  20  and  15  percent  which 
led  to  saturations  of  71  and  54  percent 
respectively . 

191.89-20-4  Run  up  to  93-g  prior  to  testing  at  26-g 

(Balance  light  did  not  come  on) . 

There  was  a  high  range  within  the  moisture 
content  measurements.  Moisture  contents 

ranged  from  10  to  18%.  Saturations  ranged 
from  38  to  64%. 


120.63-60-1 


102.08-40-2 


191 . 89-60-3 


164.71-20-5 


Eliminated  one  point;  A  moisture  content  of 
3.31%  mid  way  down  the  specimen  was 
inconsistent  with  the  moisture  content  data 
above  and  below  it . 

Bottom  measurements  are  actually  2  to  4  inches 
off  the  bottom  of  the  centrifuge  bucket. 

More’  water  was  measured  in  the  center  of  the 
specimen  than  by  the  side  of  the  centrifuge 
bucket.  Moisture  content;  20.6%  vs.  17%; 
Saturation:  74%  vs.  60%. 

Also,  the  top  and  middle  points  of  the  12 
o'clock  sample  location  were  switched. 

Moisture  content  measurements  were  larger  at 
the  center  of  the  specimen  than  along  the  side 
of  the  centrifuge  bucket.  The  top 
measurements  showed  moisture  contents  of  4.6% 
at  the  center  and  2.65%  at  the  side  which 
correspond  to  saturations  of  16.4%  and  9.5% 
respectively . 

Moisture  content  measurements  were  larger  at 
the  center  of  the  specimen  than  along  the  side 
of  the  centrifuge  bucket.  The  bottom 
measurements  showed  moisture  contents  of  12.2% 
at  the  center  and  7.6%  at  the  side  which 
correspond  to  saturations  of  43.8%  and  27.4% 
respectively . 
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TABLE  D1  SPECIMEN  MOISTURE  CONTENT  DATA 

TYNDALL  BEACH  SAND  -  DRY  DENSITY  OF  95tVfl**3  (1521.7kg/m*‘3) 


2  to 


<X)  5 

^  g  CO  n  ^  j;]  |cy  Icj  oi  CO  CM  IT)  ^  jco  ^  CO  |(M  iCM  O)  CO 


*“  CM  CM  CO 


;S^ci<i66S6 

wcocooooooo'ssoo 

CD  '  O  O  O  O  >  O 
^  Ecsic\ic\jcsi.^CNi 
9-0  0  o  o  *  o 


Attemale  depuhs  measured  from  Ihe  top  of  the  cenfrifix,^  bucket. 


TABLE  D.2  SPECIMEN  MOISTURE  CONTENT  STATISTICS 
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Regression  statistics  were  not  run  on  test  data  having  five  sample  depths  of  1  -2",  3-4”,  5-6”,  7-8",  and  9-10". 
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WATER  RETENTION  (DESATURATION)  CURVE 
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Capillory  Pressure  (cm) 


Figure  E.l  Water  Retention  ( Decaturat ion )  Curve  for 
Tyndall  Beach  Sand. 
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Copilfory  Pressure  (KPa) 


APPENDIX  F 

ACTUAL  DRY  DENSITY  OF  PLUVIATED  DRY  SPECIMENS 
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TABLE  F.1  ACTUAL  DRY  DENSITIES  FOR  THE 
PLUVIATED  SPECIMENS. 


Test  no. 

Measured 

Weight 

Bucket 

Weight 

Bucket 

Volume 

Actual 

Density 

Actual 

Density 

(lb) 

(lb) 

m) 

(lb/ft^3) 

(kq/m''3) 

1 

196.13 

29.10 

1.76 

94  96 

1521.04 

2 

197.38 

29.75 

1.77 

'34.65 

1516.15 

3 

195.75 

29.10 

1.76 

94.74 

1517.63 

4 

195.50 

29.75 

\T7 

93.59 

1499.19 

5 

195.75 

29.10 

1.76 

94.74 

1517.63 

6 

196.50 

29.75 

1.77 

94.16 

1 508  23 

7 

198.63 

29.75 

\T7 

95.36 

1527.45 

8 

195.00 

29.10 

1.76 

94.32 

1510.80 

9 

198.38 

29.75 

1.77 

'35.21 

1525.19 

10 

195.50 

29.10 

1.76 

94.60 

1515.35 

11 

195.38 

29.10 

1.76 

94.53 

1514.21 

Only  first  eight  tests  reported. 
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APPENDIX  G 


APPARENT  CRATER  MEASUREMENT  PROCEDURE 


APPARENT  CRATER  MEASUREMENT  PROCEDURE 


A.  Inside  Diameter,  and  Depth 

The  apparent  crater  depth  and  inside  diameter  were 
measured  immediately  following  the  test.  All  measurements 
were  taKen  using  the  original  ground  surface  as  a  frame  of 
reference.  A  Wang  word  processing  ruler  having  metric  units 
was  placed  across  "the  crater  and  pushed  down  through  the 
crater  lip  to  the  original  ground  surface.  The  apparent 
crater  diameter  was  measured  as  the  distance  along  the  ruler 
bounded  by  the  intersections  of  the  ruler  with  the  crater  wall 
(Figure  G.l).  Using  a  second  ruler,  the  apparent  crater  depth 
was  measured  as  the  distance  from  the  original  ground  surface 
as  established  by  the  Wang  ruler  to  the  deepest  point  of  the 
crater  (Figure  G.l). 

Indentations  left  by  the  Wang  ruler  are  shown  in  Figure 
G.2.  The  two  indentations  shown  in  the  North-South  direction 
indicate  the  distance  between  the  observed  center  of  the 
apparent  crater  and  the  location  of  the  deepest  point  in  the 
crater.  For  the  crater  shown,  the  diameter  measurement  in  the 
North-South  direction  was  made  at  the  crater  center  along  the 
eastern  most  North-South  indentation,  while  the  depth 
measurement  was  made  along  the  western  most  North-South 
indentation . 

The  apparent  crater  depth  and  diameter  measurements  were 
taken  in  two  perpendicular  directions.  The  directions  of 
measurement  were  those  providing  the  greatest  accuracy.  No 
attempt  was  made  to  limit  crater  measurement  to  specific 
directions  or  orientations.  Thus,  the  dimensions  reported  in 
Appendix  B  are  an  average  cl  the  two  measurements  taken  for 
each  crater. 

B .  Volume 

Previously,  crater  volume  measurements  have  been 
made  using  surveying  instrumentation  or  a  prof ilometer .  Both 
techniques  measure  distances  from  an  arbitrary  datum,  set 
parallel  to  the  original  ground  surface,  to  different 
elevations  within  the  crater  and  approximate  crater  volume 
using  concentric  disks.  The  result  is  a  topographic 
representation  of  the  crater. 

Although  these  techniques  are  sufficiently  accurate, 
an  alternative  method  of  measuring  crater  volume  was  designed. 
A  concrete  mold  was  made  of  the  crater  which  was  then  used  to 
make  a  clay  replica.  When  the  clay  replica  was  filled  with 
water,  the  weight  of  water  in  grams  equaled  the  crater  volume 
in  cubic  centimeters  (Figure  G.3). 

Permoplast,  artist  modeling  clay,  was  used  to  make  the 
clay  replicas  (Figure  G.4) .  Saran  wrap  was  used  to  keep  the 
clay  from  sticking  to  the  concrete  mold.  After  being  molded 
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Wang  Ruler 


Figure  G.l  Apparent  Crater  Depth  and  Inside  Diameter 
Measured  using  the  Original  Ground  Surface 
as  a  Frame  of  Reference. 
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North 
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Figure  G.3  Clay  Replica  Filled  With  Water. 

Apparent  Crater  Volume  Assumed  to 
be  the  Weight  of  Water 
{1  gram  water  =  1  cubic  centimeter). 


to  the  shape  of  the  apparent  crater,  the  clay  was  covered  with 
masking  tape  to  prevent  the  shape  of  the  crater  from  changing 
when  the  concrete  mold  was  removed  (Figure  G.5). 

Rarely  could  the  apparent  crater  depth  and  inside 
diameter  be  duplicated  when  the  clay  replica  was  filled  with 
water  even  when  the  indentations  left  by  the  Wang  ruler  were 
used  as  reference  points.  Either  one  of  the  apparent  crater 
dimensions,  diameter  or  depth, could  be  matched,  but  matching 
both  dimensions  at  once  often  proved  impossible.  Therefore, 
two  apparent  crater  volume  measurements  were  made  for  every 
concrete  mold. 

The  first  measurement  was  taken  at  the  apparent  crater 
dimension  perceived  as  being  the  most  accurate,  and  the  second 
measurement  was  taken  at  the  least  accurate  dimension.  Thus, 
the  second  apparent  crater  volume  measurement  established  an 
upper  or  lower  bound  for  the  first  measurement.  The 
difference  between  the  first  and  second  measurements  is  the 
potential  error  associated  with  the  apparent  crater  volume. 
The  potential  error  is  assumed  to  be  symmetric  so  that  the 
difference  between  the  first  and  second  apparent  crater  volume 
measurements  is  interpreted  as  a  95  percent  confidence 
interval . 


Figure  G.5  Clay  Replica  Covered  With  Masking  Tape 
to  Prevent  the  Shape  of  the  Crater  From 
Changing  When  the  Concrete  Mold  was 
Removed . 
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CRATER  MOLD  MEASUREMENT  PROCEDURE 


1.  The  crater  is  leveled  such  that  three  reference 
points  all  have  approximately  the  same  vertical  distance, 
from  a  horizontal  datum  (Figure  H.l) .  These  reference  points 
are  located  on  the  edge  of  the  mold  beyond  the  crater  lip. 
They  should  be  close  to  the  original  sand  surface. 

The  purpose  for  leveling  the  crater  mold  is  to 
attain  the  same  orientation  present  when  it  was  removed  from 
the  specimen.  In  this  manner,  the  angle  0  can  be  measured. 

2.  The  ruler  is  leveled  to  insure  the  measurements  are 
indeed  vertical. 

3.  Using  the  arrow  scored  on  the  underside  of  the  mold, 
the  ruler  is  oriented  perpendicular  to  the  rotor.  This  is  to 
place  the  profile  along  the  East-West  direction,  the  line  of 
Coriolis  acceleration. 

4.  Record  the  location,  x  distance  of  each  of  the 
reference  points. 

5.  Calculate  the  center  of  the  center  of  a  symmetrical 
crater  by  dividing  the  distance  between  reference  points  in 
half. 

6.  Measure  the  vertical  distance  from  the  ruler  to  the 
crater  mold  at  intervals  between  the  reference  points. 

7.  Calculate  the  difference,  x  distance,  between  the 
location  of  the  maximum  unskewed  crater  depth  and  the  location 
of  the  maximum  skewed  crater  depth. 

8.  Determine  the  angle  of  rotation,  (!>. 

tan<|)  =  distance  separating  maximum  crater  depths 
maximum  skewed  crater  depth 
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APPENDIX  I 

APPARENT  CRATER  DIMENSIONS  VERSUS  CHARGE  MASS,  TNT 
EQUIVALENT  (FIGURES  8.1  -  8.6)  WITHOUT  1  G  DATA 
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Scaled  Explosive  Moss.  TNT  Equivalent  (kg) 
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Figure  1.3  Scaled  Apparent  Crater  Depth  Versus  Scaled  Explosive  Mass 
TNT  Equivalent,  (Conversions:  0.305ft  =  2.0m, 

0.4541b  =  1 .0kg)  . 


Scolcd  £xp1osiv«  Moss,  TNT  Equlvolcnt  (kg) 


Scaled  Apparent  Crater  Volume  Versus  Scaled  Explosive  Mass 
TNT  Equivalent.  (Conversions:  O.C283ft^  =  l.Om^, 


Scaled  Explosive  Mass.  TNT  Equivalent  (Wg) 
10  1  10  10  ’  10  ’  1C 
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Figure  1.5  Scaled  Apparent  Crater  Radius  Versus  Scaled  Explosive  Mass 
TNT  Equivalent.  (Conversions;  0.305ft  =  1.0m, 
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Figure  1.6  Scaled  Apparent 


TABLE  1.1  YIELD  EXPONENTS  FOR  BEST-FIT  LINES  WITHOUT  1  G 
DATA. 


Dry  (S=0% 


Volume 

0.831 

0.813 

Fladius 

0.297 

0.289 

Depth 

0.279 

0.159 

TABLE  1.2  CORRELATION  COEFFICIENTS  FOR  BEST-FIT  LINES 
WITHOUT  1  G  DATA. 


APPENDIX  J 


TORQUE  P 


ON  THE  CRATER  SOIL  MASS  BY  THE  CORIOLIS  FORCE 


TORQUE  PLACED  ON  THE  CRATER  SOIL  MASS  BY  THE  CORIOLIS  FORCE 


Figure  J.l  shows  the  torque  placed  on  the  ejected  soil 
mass  by  the  Coriolis  force  (F^) .  Given  a  spherical  control 
volume  of  radius  R  centered  at  the  detonator,  the  torque  (T) 
at  any  point  can  be  expressed  as  the  triple  integral: 

f  =  J j j fR  X  df  d6  f 9in4>  ci|)  (J.l) 


where  rR  is  the  particle  velocity,  0  is  the  angle  between  r 
and  the  z  axis,  and  6  is  the  angle  between  the  trace  of  r  in 
the  x-y  plane  and  the  x  axis.  The  Coriolis  force  is 

( J.  2) 


which  can  be  rewritten  as 


F^  =  X  Vp(f) 


(J.3) 


or 


Fg  =  X  Vpdcosasintj)  +  jsinesin4) 


.^cos<|>) 


(J.4 

) 


where  m  is  a  unit  mass,  a^  is  the  Coriolis  acceleration,  O)  is 
the  angular  velocity,  Vp  is  the  particle  velocity,  and  i,  j, 
and  k  are  unit  vectors  in  the  x,  y,  and  z  directions 
respectively.  Thus,  the  torque  at  radius  R  becomes 


71  2ll 

r  =  j*  J'zimiiVplJc  (coBdsin<t>}^  +  ^(sin0sin<j>)  ^ 
0  0 


+  ( -icos4)coB0sin(J> 

-  jcos4»Gin0sin4>)  3  d0  fsin4)  di()  (j  5) 

The  integration  of  the  torque  with  respect  to  9  from  0  to 
2rt  produces  a  rotational  vector  for  which  the  X  and  Y 
components  cancel  out.  By  definition,  the  torque  is 

perpendicular  to  and  F^.,  so  the  plane  containing  T  and  Vp 
passes  through  the  origin.  For  every  velocity  vector,  Vp, 
there  is  a  second  velocity  vector  with  an  equal  and  opposite 
X  and  Y  components. 
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APPENDIX  K 

CENTRIFUGE  CRATER  ROTATION  VERSUS  G-LEVEL 
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APPENDIX  L 

DECLINE  OF  NORMALIZED  APPARENT  CRATER  DIMENSIONS 
WITH  INCREASED  G-LEVEL,  SCHMIDT  AND  HOLSAPPLE 
(1980),  AND  GOODINGS  ET  AL.  (1988). 
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An  analysis  of  the  normalized  apparent  crater  depth  data 
presented  in  Figure  8.11  was  made  to  determine  the  separation 
observed  in  the  data.  The  results  of  the  analysis  are  shown 
in  Table  L.l.  The  normalized  apparent  crater  depth  decreases 
with  an  increase  in  G-level.  Figures  L.l  through  L.4  add  the 
data  reported  by  Schmidt  and  Holsapple  (1978)  and  Goodings  et 
al.,  (1988)  to  Figures  8,9  and  8.10. 


Table  L.l  NORMALIZED  APPARENT  CRATER  DEPTH  DATA. 


Apparent  Crater  Depth/W''l/3 
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(m/kg''l/3) 
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0.80 

26 

0.81 

49 

0.60 

67 

0.57 

35,  40 

19 

0.77 

26 

0.94 

49 

0.53 

67 

0.45 

53,  60 

19 
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26 

0.75 

49 

0.53 

67 

0.51 

70 

19 

0.78,  0.78 

26 
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49 
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Figure  L.l  Normalized  Apparent  Crater  Diameter  Data 
Reported  by  Schmidt  and  Holsapple  (1978) 
Plotted  in  Figure  8.9. 
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Figure  L.2  Normalized  Apparent  Crater  Depth  Data 

Reported  by  Schmidt  and  Holsapple  (1978) 
Plotted  in  Figure  8.10. 
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Figure  L.3  Normalized  Apparent  Crater  Dianete 
Reported  By  Goodings  et  al .  (1983} 

Plotted  in  Figure  8.9. 
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